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ABSTRACT 


An  equation  of  state  model,  including  the  solid  to  liquid  and 
liquid  to  vapor  phase  transitions,  has  been  developed.  It  is 
applicable  for  materials  which  usually  melt,  rather  than  subli¬ 
mate,  and  whose  vapor  consists  of  atomic  species.  The  parameters 
appearing  in  the  formulation  are  completely  determined  from  a 
limited  amount  of  data  which,  for  many  materials,  is  currently 
available.  A  computer  program  which  gives  the  pressure  as  a 
function  of  volume  and  internal  energy  has  been  written  and  used 
to  obtain  numerical  results  for  aluminum,  titanium,  and  beryllium. 
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SECTION  1 

INTRODUCTION 


The  problem  of  stress  wave  generation  and  propagation  resulting  from  the 
deposition  of  a  relatively  large  amount  of  energy  in  a  very  short  time 
has  been  of  considerable  interest  for  over  fourteen  years.  Various 
numerical  techniques  have  been  developed  for  solving  the  coupled  set  of 
hydrodynamic  equations  with  a  high-speed  digital  computer*  Numerous 
versions  of  a  cme -dimensional  hydrodynamic  code,  generally  referred  to 
as  PUFF,  have  evolved  and  there  has  been  a  continuing  effort  to  improve 
their  capability  for  accurately  predicting  the  properties  of  the  ensuing 
stress  wave  and  its  effect  on  various  mater  in  is*  One  objective  of  the 
UASA  sponsored  PREDIX  Program  (formerly  referred  to  as  the  SUPERPUFF 
Program)  has  been  to  provide  better  input  information  for  such  a  computer 
code . 

While  the  basic  one -dimensional  hydrodynamic  equations  are  relatively 
simple,  the  accuracy  of  the  results  of  any  numerical  computation  is  no 
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better  than  that  with  which  one  can  describe  the  eneruv  deposition  process 
and  the  equation  of  state  for  the  material  of  interest.  Xhi3  report 
summarizes  Lhe  results  of  an  investigation  pertaining  to  the  equation  of 
state  problem. 

The  simplest  equation  cf  state  which  describes  the  solid  phase  for  small 

strains  and  low  energy  densities  is  the  familiar  Hooke's  Law  with  an 

added  thermal  expansion  term.  A  more  detailed  equation  is  one  which 

fares  a  modified  form  of  this  equation  into  the  perfect  gas  equation  at 

low  mass  densities  and  high  energy  densities.  While  such  an  equation 

covers  the  full  range  of  densities  and  energies  it  masks  all  details 

of  the  bolid-liquid  and  liquid-vapor  phase  transitions.  More  recently* 

an  equation  of  state  for  aluminum,  in  tabular  form,  has  been  developed 

1  2 

by  Shock  Hydrodynamics ,  Inc.  J  This  equation  of  state  is  a  modification 
of  the  McCloskey  equation  of  state  in  which  the  liquid-vapor  phase  transi¬ 
tion  is  treated  in  some  detail.  While  this  feature  is  retained  in  the 
present  formulation  the  empirical  procedures  for  describing  the  liquid 
phase  are  quite  different. 

The  equation  of  state  as  visualized  for  this  study  i'  depicted  in  Figures 
1.1a  and  1.1b  for  a  material  which  expands  on  melting.  It  is  convenient 
to  think  in  terms  of  the  thermodynamic  variables  of  pressure,  volume, 
and  temperature.  Eventually  the  internal  energy  will  be  introduced  to 
replace  the  temperature. 
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FIGURE  1.1b  P-V  CURVES  (DASHED  LINES)  FOR  THE 
MATERIAL  OF  FIGURE  1.1a 


FIGURES  1.1  EQUATION  OF  STATE  FOR  A  MATERIAL 
WHICH  EXPANDS  UPON  MELTING 


A  P-T  phase  diagram  showing  the  solid,  liquid,  and  vapor  phases  separated 
by  lines  along  which  phase  transitions  occur  is  shown  in  Figure  1.1a. 

Point  C  is  the  critical  point  represented  by  the  critical  pressure  P  , 
critical  volume  ,  and  critical  temperature  ,  while  point  A  is  the 
triple  point  where  all  three  phases  can  exist  in  equilibrium.  The  phase 
transition  lines  become  mixed  phase  regions  in  the  P-V  diagram  of  Figure 
1.1b.  One  might  expect  the  solid-liquid  region  to  have  a  top,  as  does 
the  liquid-vapor  region,  but  none  has  ever  been  observed. 

An  equation  of  state  as  shov»n  in  Figure  1.1  is  an  equilibrium  situation 
and  hence  equilibrium  is  assumed  throughout.  Rate  effects  can  be 
included  in  the  equation  of  state  for  the  solid  phase  without  affecting 
the  rest  of  the  development.  Furthermore,  the  pressure  of  Figure  1.1 
is  a  hydrostatic  pressure,  equal  in  all  directions.  For  one-dimensional 
stress  wave  propagation  in  materials  which  can  support  shears  the  stress 
normal  to  the  direction  of  propagation  is  different  from  that  in  the 
direction  of  propagation.  This  problem,  which  arises  only  in  the  solid 
phase,  will  be  considered  at  the  appropriate  time. 

The  equation  of  state  depicted  in  Figure  1. 1  does  not  allow  for  negative 
pressures  in  the  solid  phase,  and  the  solid  sublimates  at  sufficiently 
low  pressure.  For  the  materials  considered,  however,  the  available  data 
indicates  that  the  triple  point  pressure,  if  a  triple  point  exists,  is 
many  orders  of  magnitude  below  one  atmosphere.  Hence  it  is  assumed  that 
the  pressure  at  point  A  can  be  taken  as  zero  and  therefore  the  sublimation 
region  is  neglected.  Furthermore,  the  melting  line  AB  is  extremely  steep 


so  the  temperature  at  point  A  can  be  taken  eaual  to  the  normal  melting 
temperature.  The  resulting  modified  equation  of  state  is  shown  in 
Figure  1.2.  Negative  pressures  in  the  solid  phase  are  now  possible  for 
temperatures  less  than  the  normal  melting  temperature  or  for  internal 
energies  less  than  that  at  which  the  material  begins  to  melt.  This 
simplification  does  not  apply,  of  course,  to  materials  like  carbon  which 
normally  sublimate. 

The  problem  then  is  to  determine  parameters  for  the  equation  of  state  of 
Figure  1.2  for  materials  of  interest.  Unfortunately,  the  amount  of  data 
for  most  materials  is  quite  limited  and  one  is  forced  to  make  numerous 
assumptions.  An  attempt  is  made,  however,  to  keep  the  formulation  some¬ 
what  flexible  so  that  any  additional  data  which  might  become  available 
can  be  included  to  refine  the  parameters. 

Section  2  summarizes  the  available  data,  not  only  for  the  materials 
considered  in  the  study,  but  also  for  other  materials  where  experimental 
measurements  around  the  critical  point  have  been  made.  This  enables  one 
to  retain  certain  qualitative  features  of  an  equation  of  state  where 
data  is  known  but  also  indicates  marked  differences  between  different 
types  of  materials. 

The  bulk  of  the  effort  is  contained  in  Section  3  where  analytical  expressions 
for  various  regions  of  the  PV  diagram  are  developed.  The  liquid  phase  is 
the  mo.it  difficult  to  describe  and  the  procedure  taken  is  to  approach  it 
both  from  the  vapor  phase  through  the  vapor-liquid  transition  region 
and  from  the  solid  phase  through  the  solid-liquid  transition  region.  In 
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PRESSURE 


FIGURE  1.2a  MODTFIF.il  PHASE  DIAGRAM 


VOLUME 


FIGURE  1.2b  MODIFIED  P-V  CURVES  (DASHED  LINES) 
FOR  THE  MATERIAL  OF  FIGURE  1.2a 


FIGURE  1.2  MODIFIED  EQUATION  OF  STATE  FOR  A  MATERIAL 
WHICH  EXPANDS  UPON  MELTING 
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the  process  of  developing  the  analytical  expressions,  a  possible  procedure 
for  determining  the  necessary  from  a  limited  amount  of  data 

is  outlined. 

Section  4  describes  an  equation  of  state  subroutine  which,  given  a  volume 
V  and  an  energy  E,  determines  the  pressure  P.  An  initial  calculation 
determines  the  required  parameters  from  relatively  simple  input  data. 

The  results  of  some  calculations  using  this  subroutine  are  given  in 
Section  5  for  aluminum,  titanium  and  beryllium.  Tsoenergy  curves  and 
adiabats  have  been  determined  to  check  out  the  subroutine. 
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SECTION  2 

DATA 


The  amount  of  data  for  various  materials  is  quite  extensive  for  those 
which  are  vapors  under  normal  conditions  or  become  vapors  at  very  low 
temperatures,  rather  complete  for  materials  such  as  mercury  and  the  alkali 
metals  which  vaporize  at  reasonable  temperatures,  and  very  limited  for 
metals  whose  boiling  temperatures  are  quite  high.  This  section  sum¬ 
marizes  some  of  the  data  which  is  useful  as  a  guideline  in  the  subsequent 
development.  The  more  pertinent  data  is  plotted  in  Figures  2.1  and  2.2. 

2.1  MATERIALS  WITH  LOW  BOILING  TEMPERATURES 

The  data  for  materials  which  are  vapors  under  normal  conditions  and  those 

which  boil  at  very  low  temperatures  have  been  analyzed  in  detail  by 

3  4  5 

Hirschf e lder ,  et.  al.  9  9  These  include  the  noble  gases,  hydrocarbons , 
and  the  highly  polar  substances,  water  and  ammonia.  Data  exists  for  rather 
wide  ranges  of  temperature,  density,  and  pressure  both  below  and  above  the 
critical  point  values.  The  result  of  the  analysis  is  a  generalized 
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equation  of  slate,  containing  relatively  few  parameters,  lor  the  vapor, 
liquid,  and  mixed  phase*. 

An  important  parameter  which  appears  throughout  is  the  critical  compressi¬ 
bility  factor  Z defined  by 

Z  =  P  V  /K  T  (2.1) 

c  c  c  c 

where  P^,  ,  and  are  the  critical  pressure,  molar  volume,  and  tempera¬ 

ture,  respectively.  K  is  the  molar  gas  constant.  Typical  values  of  Z 
for  these  materials  are: 

Substance 

Noble  gases  0.29-0.30 

Light  hydrocarbons  0.27 

Water  0.23 

While  substances  which  vaporize  at  high  temperatures  do  not  fall  into 
the  same  regime  as  those  analyzed  by  Hirschfe lder ,  et.  al,  several  of 
the  analytical  forms  proposed  by  them  have  been  retained  to  describe  the 
vapor  phase  and  the  liquid -vapor  mixed  phase. 

2 . 2  MERCURY 

Data  on  the  thermodynamic  properties  of  mercury  are  more  complete  and 

probably  more  accurate  than  for  any  other  metal.  Measurements  of  the 

6  7  8 

vapor  pressure  9  and  the  saturated  liquid  and  vapor  densities  have  been 

made  to  very  near  the  critical  point.  On  the  basis  of  these  data,  the 

9 

critical  parameters  are  estimated  to  be 
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1587  atm 


T  -  J  73  J  K 

c 


6.7  gm/cc 


2  -  0.48 

c 


By  assuming  a  principle  of  corresponding  states,  that  the  vaporization 

9 

entropy  of  various  liquids  *a  equal  at:  corresponding  temperatures,  Crosse 
also  estimates  the  critical  temperatures  for  other  materials.  This  idea 
is  retained,  but  some  weight  is  a 'so  given  to  the  data  for  the  alkali 
meta  Is . 

2. '3  ALKALI  METALS 

The  vapor  pressures  of  the  alkali  metals  have  been  measured  up  to  100  psi^ 
and  the  data  this  lov;  temperature,  low  pressure  range  can  be  fit  with  a 
Raou 1 L  relation  of  the  form 

log  P  -  A  -  B/T.  (2.2) 

The  constants  A  and  B  are  given  in  Table  2.1,  when  the  temperature  is  In 
degrees  Kelvin  and  the  pressure  is  in  atmospheres. 

Table  2.1:  Vapor  Pressure  Constants  for  the  Alkali  Metals 

Element  A  B 


4.52172 


5220.42 


4.09537 


4205.78 


4.04261 


3880.020 


3. 941504 


3702.814 


Dillon,  et.  al,  have  measured  the  saturation  densities  of  both  liquid 
and  vapor  up  to  near  the  critical  point  and  from  this  data  they  estimate 
the  critical  properties  of  the  alkali  metals.  The  critical  parameters 
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are  listed  in  Table  2.2. 


Table  2.2:  Critical  Parameters  tor  the  Alkali  Metals. 

Atomic  Boiling  Critical  Critical  Critical  Coroprp««i M  H t v 

Element  WeiRht  Temp . (  K)  Temp.  1  K)  Density  (um/cc)  Press .  (atm)  Fact or 

Na  22.99  1183  2573+150  0.206  +0.041  350  0.20+0.12 

K  39.10  1039  2223+330  0.194  +0.037  160  0.21+0.13 

Kb  85.47  974  2093  +  35  0.346  +  0.009  157  0.22  +  0.02 

Cs  132.90  958  2057  +  40  0.428  +  0.012  145  0.25  +  0.02 

^Estimated  Error  +  207., 


The  values  shown  were  taken  from  their  table  and  there  is  no  explanation 
for  slight  discrepancies  which  appear.  In  the  following  it  is  assumed  that 
the  critical  temperatures  and  densities  are  correct. 

This  data  can  now  be  used  to  construct  curves  of  vaporization  entropy  as 
a  function  of  the  reduced  temperature,  t  *  T/T^,  at  least  for  the  lower 
temperatures,  in  the  following  manner.  The  Clapeyron  equation  for  the 
slope  of  the  vapor  pressure  curve  is 


(2.3) 


where  /.H  is  the  latent  molar  heat  of  vaporization  at  the  temperature  T  and 


pressure  P,  Vv  is  the  molar  volume  of  the  saturated  vapor,  and  is  the 

molar  volume  of  the  saturated  liquid.  At  low  temperatures,  ,  and 

it  is  further  assumed  that  V  is  that  for  an  ideal  gas.  Hence 

v 


V  -  V.  V  RT/P. 
v  i  ‘  v  - 


(2.4) 


2  ”4 


The  vaporization  entropy  i«  then 


fil  df 

T  ^  P  dT  ' 


iu  U.3U/b 


“)  -  2  .  3026 


=  2.3026 


where  use  has  been  made  of  Fq .  (2.2)  to  evaluate  the  derivative.  This 

function  is  plotted  in  Figure  2.1,  using  the  value3  of  IJ  and  T  from 

Tables  2.1  and  2.2.  Also  shown  in  Figure  2.1  are  experimental  values  for 
9  9 

Hg  and  a  range  of  values  for  the  low  boiling  temperature  substances, 
This  figure  will  be  used  to  estimate  the  critical  temperatures  for  other 


meta  Is . 


Another  quantity  of  interest  is  the  reduced  saturated  liquid  density  as 

a  function  of  the  reduced  temperature.  This  is  plotted  in  Figure  2.2,  In 

addition  to  the  alkali  metal  data,  some  Hg  data  points  and  a  range  of 

3 

values  for  hydrocarbons  are  also  shown.  Hirschfelder  suggests  that  the 
data  for  a  given  material  can  be  fit  with  the  empirical  form 

1  /3 

-  1  +  c  (l-t)1/J  +  d  (1-t),  (2.6) 

where  is  the  reduced  density  and  t  is  the  reduced  temperature.  The 
alkali  metal  data  can  all  be  fit  with  the  same  curve  with  c  *  1.8  and  d  =  1,75. 
This  function  is  also  shown  in  Figure  2*2 . 

2.4  ALUMINUM,  TITANIUM  AND  BERYLLIUM 

There  is  no  direct  data  available  at  high  temperatures  and  pressures  for 
aluminum,  titanium  and  beryllium.  This  information  must  be  derived  from 
other  data  which  is  at  hand.  Table  2.3  indicates  the  type  of  information 

which  is  usually  available  for  the  metals  and  lists  values  for  the  three 
metals  being  considered. 
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I  CURE  2.1  VAPORIZATION  ENTROPY  VERSUS  REDUCED  TEMPERATURE 


REDUCED  LIQUID  DENSITY 


V 


ALKALI  METAL 
FITTED  CURVE : 


4 


V 

.  3 


2 


1 


(j 

REDUCED  TEMPERATURE,  t 


FIGURE  2.L*  REDUCED  LIQUID  DKNSTi'Y  VERSUS  REDUCED 
TEMPERATURE 
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Table  2.3.  Data  lor  Aluminum,  Titanium  a  ml  lie  ry  Ilium 


A  1 

*T  4  ~~  J 

i  um 

r» -  — i  »  J  — 
UV.  L  jr  1  11  lif.l 

Atomic  Weight  (gm/mole) 

26.98 

47.90 

9.013 

Solid : 

Normal  Density  (gm/cc) 

2.71 

4.5 

1.85 

C  (kbar)*~ 

772.0 

1016 

1203 

12 

D  (kbar) 

490.8 

"’22.2 

821.2 

S  (kbar 3^ 

607.6 

-568.5 

-379.0 

12 

C  (kbar /ca 1/gm) 

0.2384 

0.2044 

0.08938 

Poisson's  Ratio 

1/3 

Solid  to  Liquid: 

Tielting  Temperature  (°K)  ^ 

932 

1950 

1556 

13 

Solid  Enthalpy  (kcal/mole) 

4.291 

13.390 

7.910 

13 

Liquid  Enthalpy  (kcal/mole) 

6.838 

17.094 

10.714 

13 

Vapor  Enthalpy  (kcal/mole) 

81.177 

121.274 

84.504 

Heat  of  Fusion  (kcal/mole) 

2.547 

3.704 

2.798 

Solid  Density  (gm/cc) 

2.53714 

(4.25) 

1.808 15 

Liquid  Density  (gm/cc) 

2.38014 

(3.97; 

1.69015 

Slope  of  FT  Curve  (atm/°K) 

161 

f  </:  .  ?  ^ 

213 

Liquid  to  Vapor: 

o  13 

Boiling  Temperature  (  K) 

13 

Liquid  Enthalpy  (kcal/mole) 

13 

Vapor  Enthalpy  (kcal/mole) 


2736 

3550 

2757 

19.466 

29.894 

19.342 

90. 149 

132.348 

90.477 

70.683 

102.454 

71.  135 

1.92  (2 

.07) 

(3.65) 

1.551 

25.83 


(1.73) 

(1.62) 

(210) 


Heat  of  Vaporization  (kcal/mole) 
Liquid  Density  (gm/cc) 

Vaporization  Entropy  (cal/mole/  K) 
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28.86 


25.80 


(1.47) 


It  is  assumed  that:  an  analytical  form  of  the  equation  of  state  for  the  solid 
exists ,  i.e.,  pressure  as  a  function  of  density  P  and  internal  energy  E. 
Various  forms  have  been  proposed  but,  for  the  purpose  of  this  discussion, 
it  i  s  take  *i  lo  be^ 

P(P,E)  =  Cn  +  Du2  ■»  S^3  +  G  TJ  E.  (2.7) 

7?  *  P/Pq,  °0  being  the  normal  density,  and  m-  *  V  -  1.  The  internal  energy 

is  referenced  relative  to  nomal  density  and  room  temperature.  The  values 

12 

for  the  coefficients  were  taken  from  the  RADS  report.  Analytical  forms 
different  from  that  given  in  Eq.  (2.7)  can  be  handled  with  no  more  difficulty 
by  making  minor  changes  in  the  computer  program  described  later. 

For  many  metals  there  is  a  solid  phase  change  and  this  should  be  included 
in  the  solid  equation  of  state.  In  this  case  it  is  necessary  to  have  an 
expression  for  the  material  at  high  internal  energies,  around  the  energy  at 
which  it  begins  to  melt,  in  order  to  describe  the  melting  process. 

The  solid- to- liquid  and  liquid-to- vapor  phase  transitions  are  described  by 
values  of  certain  thermodynamic  quantities  at  the  normal  melting  temperature 
and  normal  boiling  temperature,  respectively.  In  particular,  it  is  assumed 
that  the  values  of  the  enthalpy  at  these  temperatures,  as  indicated  in  the 
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This  is  just  the  information  that  one  has  frequently  used  to  arrive  at  a 

total  ynertfv  nf  vi»rw>r  i  vat  i  nn  .  Th«»  #»t>hhja  lny  wa  {nos  listed  i  r!  Tsb  ifi  2*3 

13 

were  taken  from  the  JANAF  Tables* 

In  some  cases  values  of  the  solid  and  liquid  densities  at  the  melting 
temperature  have  been  measured*  If  so,  they  can  be  used  as  input  data. 

There  is,  however,  a  consistency  condition  that  the  solid  equation  of 
state  evaluated  at  the  solid  density  and  the  energy  at  which  the  solid 
begins  to  melt  yield  zero  pressure. 

In  the  case  of  aluminum,  the  solid  density  has  been  assumed  to  be  correct 
requiring  the  internal  energy  at  the  melting  temperature  to  be  increased 
by  1.44  kcal/mole  (53.52  cal/gm).  All  other  internal  energies  are  in¬ 
creased  accordingly,  retaining  the  indicated  heats  of  fusion  and  vapori¬ 
zation  which  a  -  taken  to  be  correct. 

For  titanium,  these  densities  have  not  been  measured  and  hence  it  has 
been  assumed  that  the  solid  internal  energy  at  the  melting  temperature 
(equal  to  the  enthalpy  since  the  pressure  is  taken  to  be  zero)  is  correct. 
The  solid  density  is  then  found  to  be  4.25  gm/cc,  given  in  parenthesis. 

The  liquid  density  has  been  taken  equal  to  6.57,  less  than  the  solid  density. 
This  is  consistent  with  the  density  decrease  for  other  materials  where 
data  is  available. 

Beryllium  is  unique  in  that  at  a  temperature  of  1250°C  the  crystal  structure 
changes  from  a  hep  to  a  bee  form  with  a  corresponding  increase  in  density 
of  4.87?.^  The  hep  form  has  a  density  of  1.725  gm/cc  while  the  bcc  form 
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has  a  density  of  1.808.  This  accounts  for  the  relatively  high  solid 
density  given  in  the  table.  If  this  density  is  used,  the  solid  internal 
energy  must  be  decreased  by  5.142  kcal/mole  (570.5  cal/gm).  The  value 
of  1.73  gm/cc,  given  in  parenthesis,  is  obtained  by  assuming  the  solid 
internal  energy  is  correct.  The  liquid  density  is  again  taken  to  be 
6. 57  less  than  the  solid  density.  There  are  discrepancies  of  about  4.37 
in  the  two  sets  of  values. 

The  slope  of  the  P-T  melting  curve  can  be  evaluated  from  the  Clapeyron 
equation, 

m  _ Hi _  (2  ov 

dT  T(V£  -  Vg)  9  K  ' 

where  AH,  the  heat  of  fusion,  is  equal  to  the  difference  in  the  liquid 
and  solid  enthalpies.  Calculated  values  for  this  slope  are  given  in 
Table  2.3.  Note  that  they  are  extremely  steep  (~  100-200  atm/°K) .  Hence, 
going  from  one  atmosphere  pressure  to  zero  pressure  changes  the  melting 
temperature  by  less  than  a  hundredth  of  a  degree  and  the  melting  tempera¬ 
ture  values  given  are  certainly  valid  at  zero  pressure. 

The  liquid  density  at  the  boiling  temperature  is  known  for  some  materials, 
or  can  be  estimated  by  extrapolating  liquid  density  data  above  the  melting 
temperature.  This  information  has  not  been  used,  however,  and  the  values 
given  in  parenthesis  are  those  obtained  by  the  procedure  to  be  outlined 
later.  In  the  case  of  aluminum  the  discrepancy  is  about  87„,  but  then  the 
density  value  of  1.92  gm/cc  is  highly  questionable.  It  is  less  for 
beryllium. 
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Values  for  the  vaporization  entropy,  equal  to  the  heat  of  vaporization 
divided  by  the  boiling  temperature,  are  also  given. 

At  this  point  the  critical  temperatures  can  be  estimated.  It  is  assumed 
that,  for  the  metals,  the  vaporization  entropies  at  the  normal  boiling 
temperatures  fall  in  the  range  defined  by  the  alkali  metals.  Hence, 
from  Figure  2.1,  reasonable  values  of  the  critical  temperatures  are 
8000°K,  11,600°K,  and  8000°K  for  aluminum,  titanium,  and  beryllium,  res¬ 
pectively.  This  procedure  effectively  says  that  the  critical  temperature 
for  a  metal  is  proportional  to  its  heat  of  vaporization  at  the  boiling 
temperature  and  an  approximate  formula  is 

Tc  (°K)  ...  113  Ml  (kcal/mole).  (2.9) 

If  one  side  or  the  other  of  the  alkali  range  is  chosen,  instead  of  the 
center,  the  critical  temperatures  can  vary  by  +  500°K  from  the  above 
va  lues . 

The  lack  of  high  temperature  data  makes  it  difficult  to  check  the 
assumptions  made  in  developing  an  equation  of  state  formulation.  Appendix  A 
describes  a  laser  deposition,  inertial  confinement,  experiment,  suggested 
by  the  qualitative  results  of  this  study,  which  might  provide  some  of 
this  high  temperature  information. 


SECTION  3 

EQUATION  OF  STATE 


This  section  suggests  various  analytical  forms  for  the  thermodynamic 
variables  in  different  regions  of  the  phase  diagram.  As  has  already  been 
mentioned,  the  liquid  region  seems  to  be  the  most  difficult  to  treat. 

The  procedure  adopted  is  to  approach  it  both  from  the  vapor  phase  at  very 
low  densities  and  from  the  solid  phase  at  nearly  normal  densities.  In 
view  of  the  lack  of  data,  it  is  necessary  to  make  numerous  assumptions. 

The  formulation  which  is  developed  should  be  refined  to  incorporate 
any  new  data. 

2.1  VAPOR  PHASE 

3  4  5 

The  general  formulation  of  Hi rschf e lder ,  et.  al.,  *  *  for  describing 
the  vapor  phases  of  the  noble  gases  and  hydrocarbons  is  retained.  In 
this  formulation  the  proposed  equation  of  state  for  the  vapor  is  a  genera¬ 
lization  of  the  van  der  Waals  equation  and  can  be  written 

(P  +  A(T)/V2  +  A'(T)/V3)  (V-B+B'/V)  =  RT .  (3.1) 
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Hence,  van  der  Waals*  constant  A,  related  to  the  energy  of  attraction 
between  molecules,  is  replaced  by  A(T )  +  A*(T)/V,  a  function  of  both 
temperature  and  volume,  and  the  constant  B,  related  to  the  excluded 
volume  occupied  by  the  molecules,  is  replaced  by  B  -  B'/v,  B  and  B '  being 
constants » 


A  connection  between  the  various  parameters  of  Eq *  (3.5)  is  obtained  by 
the  requirement  that  the  first  and  second  partial  derivaties  of  p  with 
respect  to  p  vanish  at  the  critical  point,  p  38  p  *  t  -  1.  Assuming  for 
the  moment  that  a'(l)  *  0,  this  gives 
a(l)  =  3 

b  =  (35 2  -  6p  -  l)/3  (38-1)  (3.6a) 

b'  =  (-  -  3 ) / ( 30  -1). 
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where  g  is  a  new  parameter  introduced  for  convenience  and  uniquely  related 


to  Zc  by 

Z  =  P  V  /RT  =  5  08  -  1 )  /  (.*  +  1)J.  (3. 6b) 

c  c  c  c 

g  =  3 ,  Z =  0.375,  a  -  3,  b  -  1/3,  b*  =  0,  with  a*  -  0,  corresponds  to 
the  original  van  der  Waals  equation. 


It  should  be 

Z  of  0.3790 
c 

in  using  the 
are  given  in 


noted  here  that  Eq .  (3.6b)  indicates  a  maximum  value  for 
occurring  when  3  =*  2.535.  Hence  there  is  this  restriction 
form  of  Eq .  (3.1).  Typical  values  of  3  as  a  function  of  Z 
Table  3.1. 


Table  3.1:  Parameter  ft  as  a  Function  of  Z  . 


z 

c 

i— 

0.23 

9.24 

0.24 

8.65 

0.25 

8.  12 

0.26 

7.62 

0.27 

7.15 

c 


Z 

_ c 

a 

0.28 

6.71 

0.29 

6.30 

0.30 

5.90 

0.31 

5.52 

0.32 

5.51 

The  temperature  dependent  forms  suggested  for  a(t)  and  a'(t)  are 

a(t)  =  k  +  (S-k  )/t  -  k  +  k,/t 

0  °  °  1  (3.7) 

a ' (t )  -  \  (1  -  k  +  28  -  or)  (t  -  1/t)  =  k,.(t  -  1/t). 

2  o  2 

k  is  an  adjustable  constant  and  Hirschf e lder ,  et.  al.  use  the  value  of 
o  J 

5.5  throughout  their  calculations.  They  point  out,  however,  that  a  value 
of  4.71  might  give  a  better  fit  to  the  noble  gasas  and  hydrocarbon  data. 
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Since  typical  Z^  values  for  the  metals  seem  to  be  around  0,31,  and  since  a 
better  value  for  kQ  is  not  known,  it  has  arbitrarily  been  taken  equal  to  0 
for  obtaining  the  results  of  Section  5*  3  is  a  parameter  discussed  in  the 

next  section  and,  for  the  present,  it  is  sufficient  to  note  that  3p/^t 
evaluated  at  the  critical  point  is  equal  to 

The  internal  energy  E  in  the  vapor  phase  can  be  derived  from  the  thermodynamic 
relationship 

(%  '  T(^)v  '  p’  °-8) 

In  terms  of  reduced  variables  this  differential  equation  becomes 


Si  ■  w 


The  internal  energy  is  now  found  by  integrating  along  an  isotherm  from  in¬ 
finite  volume  or  zero  density, 


Z  RT 
c  c 


E(v  or  p,t)  -  E  (T)  -  Z  RT  f  [t  (§£ )  -  p  j  dv 

°  c  c  •'»  \  y v  J  (3. 

f  [p  -  E  (if)  ]  %  ■  ZcRIc  [-<ko  +  2V‘>»  +  k2  p2/tJ 

0  \  ,/p  P 


(3*10) 


Eq(T)  is  the  internal  energy  of  the  completely  separated  molecules  at  the 
temperature  in  question  and  an  analytical  expression  for  it  is  required. 

For  the  present  it  is  assumed  that  the  va^or  at  zero  density  is  an  ideal 
monatomic  gas  and  that  there  are  only  translational  and  electronic  contri¬ 
butions  to  the  internal  energy.  The  contribution  due  to  translation  is 
quite  simple  and 


Et(T)  -  |  RT. 


(3.11) 


1 


There  is  also  a  contribution  from  electronic  excitation  at  the  higher 
temperatures  and  this  is  given  by 


E  (T) 


1  fei  *1  e 


k'l 


if  gi  " 


;  i 

kT 


n.  12) 


g^  is  the  statistical  weight  of  the  i ^  energy  level,  or  group  of  levels 
with  approximately  the  same  energy,  and  is  the  corresponding  energy. 
This  internal  energy  contribution  is  for  one  mole  when  N  *  6.02338  x  10^*V 
mole,  Avagadro*s  number.  Energy  conversion  factors  are  required  to  give 
specific  units  for  the  internal  energy.  A  set  of  multiplicities  and 
electronic  levels  for  aluminum  are  given  in  Table  3.2. 


Table  3.2.  Multiplicities  and  Electronic  bevels  Lor  Aluminum. 


£i 

ei(ev) 

.<°K) 

H 

•I  J  (ov) 

,:i(°K) 

2 

0 

0 

6 

4.9932 

57950 

4 

0,0139 

161 

14 

5.1226 

59452 

2 

3.1426 

36472 

12 

5.2319 

60719 

12 

3.6078 

41870 

20 

5.4054 

62734 

16 

4.0541 

47050 

18 

5.0583 

63928 

2 

4.6726 

54229 

84 

5.7191 

66374 

10 

4.8265 

56014 

LIMIT 

5.9855 

69466 

internal 

energy  for 

an  ideal 

monatomic  gas 

is  then 

E  (T)  =  E 
o  o 

+  Kt  (T) 

+  Ee(T) 

(3.13) 

3-5 


where  K  is  a  reference  energy  at  0  Ck.  This  is  eventually  chosen  1 u 
give  zero  internal  energy  for  the  solid  at  normal  density  and  room  tempera- 
lure  . 


In  principle  Kqs.  (  \ .  1 1 ) -( 3 .  1 3)  could  be  used  to  determine  the  internal 
energy  lor  an  ideal  gas.  It  is  rather  complicated  and,  in  view  of  other 
uncertainties  in  the  formulation,  it  seems  that  a  simplified  expression 
might  suffice.  There  are  various  possibilities  but  the  one  used  assumes 
that  the  enthalpy  of  an  ideal  gas  can  be  approximated  by  a  linear  function 
of  tilt*  temperature  whose  slope  is  determined  by  the  enthalpies  at  the 
melting  and  boiling  temperatures.  Hence 


H  (T)  *  H  +  C  T 
o  o  P 

where  the  effective  specific  heat  is 


S  ‘ 


Hv  (V  -  W 
tb  *  tm 


(i. 14  ) 


(3.15) 


The  expression  for  the  internal  energy  of  an  ideal  gas  is  then 


Eo(T)  -  H0(T)  -  RT  =  Eo  +  C  T  (3,16) 

where  CY  =  -  R.  From  the  enthalpy  values  of  Table  2.3,  C  values 

VP  J  V 

are  2.99,  4.93,  and  2.99  cal/mole/°K  for  aluminum,  titanium,  and 
beryllium  respectively.  One  advantage  of  this  linear  form  is  that 
Fq .  (3.10)  can  be  solved  for  t  once  E  and  are  given. 


Eqs,  (3.5),  (3.10),  and  (3.16)  now  represent  the  equation  of  state  for 
the  vapor  phase  giving  pressure  and  internal  energy  as  functions  of  the 
density  and  temperature.  They  are  u^ed  up  to  the  vapor- liquid  mixed 
phase  envelope  for  t  <  1  and  up  to  the  critical  density  for  t  >  1. 
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2  2  vAPnK-nnmn  mtvpd  phasf.. 

In  order  to  describe  the  vapor* liquid  mixed  phase  region,  it  is  necessary 
to  have  an  expression  for  the  vapor  pressure  as  a  function  of  the  tempera¬ 
ture.  Over  a  rather  limited  range  of  temperatures  the  Raoult  relation 
has  frequently  been  used  and  the  vapor  pressure  equation  is  of  the  form 

in  F  -  A  -  B/T.  (3,17) 

While  this  form  may  be  sufficiently  good,  a  somewhat  more  general  form  has 
been  adopted. 

A  starting  point  for  developing  this  more  general  form  is  the  Clapeyron 
equat i on , 


dP 

dT 


V(V  -v  J, 
v  x 


(3.18) 


where  the  subscripts  v  and  i  refer  to  the  vapor  and  liquid  phases  res¬ 
pectively.  AH  is  the  heat  of  vaporization  at  the  temperature  T.  At  the 
lower  temperatures  the  vapor  volume  is  much  larger  than  the  liquid  volume 
and  nearly  given  by  that  for  an  ideal  gas.  Hence 


V  »  V  ^  V  -v  RT/P.  (3.  19) 

v  £  v 

Furthermore  the  heat  of  vaporization  at  these  low  temperatures  is  approxi¬ 
mately  given  by 


^  AH„  +  'C 


(T 


(3.20; 


where  AC  is  the  difference  in  specific  heats 
vapor  and  liquid  and  the  subscript  B  denotes 
ture.  To  the  same  approxintat ion  used  in  the 


at  constant  pressure  for  the 
the  normal  boiling  tempera- 
previous  section  the  vapor 
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sped  tic  heat  is  given  by  Eq.  (1.15)  and  the  liquid  specific  heat  is  ^ivcn 


by 


<V  *  H, 


(T„) 


(3.21) 


Eqs.  (3.19)  and  (3.20)  make  it  possible  to  integrate  Eq .  (3.18)  with 
the  result  that 


'b  •  /  CT» 


rJ  +  -Lr £n  ir-  0,2:0 


Ideally  ,  the  pressure  at  the  normal  boiling  temperature  Tu,  is  one 


B 


atmosphere  but  t h i s  condition  is  relaxed  somewhat. 


It  is  convenient  to  express  Eq.  ( li .  22)  in  terms  of  reduced  variables  and 
then  it  can  be  written 


fn  p  *  (c< -a^)  (1-1/t)  +  oi  ^n  t, 


U.2J) 


where  &  and  a  ^  are  dimensionless  parameters  defined  by 


Or 


ahb  +  AC  (Tc  -  rB)  _  ^ 

RT  ’  al  R 


(3.24) 


Eq .  (3.23)  is  now  assumed  to  hold  throughout  the  temperature  range  t  1 
with  a  and  & ^  given  by  Eq .  (3.24).  cr  is  the  same  quantity  appearing  in 

Eq.  (3.7)  and  it  is  the  slope  of  the  reduced  vapor  pressure  curve  evalu¬ 
ated  at  the  critical  point.  Using  the  critical  temperatures  deduced  in 
Section  2.4,  values  of  o'  for  aluminum,  titanium  and  beryllium  are  3.78, 
4.08  and  3.74,  respectively. 


In  principle,  one  could  assume  a  normal  boiling  pressure  of  one  atmosphere 
and,  knowing  a  critical  temperature,  determine  the  critical  pressure  from 
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Eq.  (3.22).  There  is  the  possibility  of  considerable  error  in  doing  this, 
however,  due  to  the  rapid  variation  of  pressure  with  temperature  and  the 
uncertainty  in  the  vapor  pressure  curve  around  the  critical  point.  There¬ 
fore  an  alternative  approach  which  iirst  determines  the  less  sensitive 
parameter  Z^  has  been  used. 

Hi rschf e  lder ,  et.  al.  ,  use  a  form  for  the  reduced  vapor  pressure  suggested 
by  Riedel, 


in  p  =  in  t  -  0.0838  -  3.75)  (  7^  -  35  -  t6  +  42  £11  t).  (3.25) 

K  K  t 

2  is  the  slope  of  this  vapor  pressure  curve  at  the  critical  point  and  is 
the  only  parameter  appearing,  Hence,  if  one  used  this  form  he  would  be 
inclined  to  choose  u ^  to  satisfy  the  Clapeyron  equation  at  the  boiling 
temperature*  Eq.  (3.19)  is  valid  here  and  hence 


determi nes 
between  oe 

K 


i  12. 

p  dt 


''  H, 


<  t=t , 


RT 


G  .  In  addition  to  Eq .  (3.25),  Riedel  suggests  a 
K 


and  the  critical  compressibility  factor  Z given 


(3.26) 

re lat ionship 
by 


Zc  =  1^3.72  +  0.26  (aR  *  7  J  .  (3.27) 

While  this  procedure  for  determining  is  questionable,  since  Eqs.  (3.25) 
and  (3.27)  are  based  on  hydrocarbon  data,  it  does  give  reasonable  values 
around  0.307  for  the  three  metals  being  considered. 


The  vapor  side  of  the  vapor i zat ion  envelope  is  now  determined  by  setting 
the  equation  for  the  pressure  of  the  vapor,  Eq.  (3.5),  equal  to  the  vapor 
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pressure,  Eq.  C3.23),  and  solving  for  the  vapor  density  as  a  function 
of  t  <  1, 

Cvt /Zt.  (  1  -bo  v+b  ’  0v2)  -a\L)ov2-a(t)pvJ  =  expj^(o-ai)(l-l/t)  +  <y  j  .;n  tj  (3. 28) 

In  this  equation,  all  variables  are  reduced.  There  is  always  a  solution 
since  the  form  of  the  equation  for  the  pressure  of  the  vapor  was  chosen  such 
that  “>p/”t  =  a  at  the  critical  point.  This  ensures  that ,  for  a  given 
t  <  1,  the  vapor  pressure  is  less  than  the  first  peak  in  the  modified 
van  der  Waals  equation  for  the  pressure  of  the  vapor.  The  internal 
energy  ^(t)  along  this  envelope  is  then  found  from  Eq .  (3*10). 

The  liquid  side  of  the  vaporization  envelope  must  be  determined  independent ly 
and  the  form  suggested  by  Hirschfe Ider ,  et.  al.,  is  retained  for  doing 
this.  This  assumes  that  the  reduced  liquid  density  p  can  be  adequately 

Xj 

represented  by 

p£(t)  =  1  +  c(l-t)1/3  +  d  (l-t)  (3.29) 


where  c  and  d  are  constants  for  a  particular  material.  The  internal 
energy  along  the  liquid  side  of  the  envelope  is  determined  from  its  value 
on  the  vapor  side  and  the  Clapeyron  equation,  Eq .  (3.18),  In  terms  of 
reduced  variables, 


dp  (t)  _  _ 


dt 


Z  RT 
c  c 


H  -H, 

V  A  - 

t (v  -V  .  x 

V  t) 


Z  RT 
c  c 


E  * 

v  JL 


!  t(vv-v;)  t  j 


+  5-|.  (3.30) 


and 


Ec(t)  =  Ev(t)  -  2.c  RTc(vv-v£)t  -  -f- j 

=  Ey(t)  -  ZcRTc  (vv-vf)  p(t)  [te-wp/t  +«!-!] 


(3.31) 
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Hence,  the  internal  energy,  pressure,  and  density  are  known  on  the  liquid 
side  of  the  saturation  envelope. 

In  the  vapor-liquid  mixed  phase  region  it  is  convenient  to  introduce  the 
vapor  mass  fraction  f  and  then  volume  V,  internal  energy  E,  and  f  are 
related  through  the  temperature  t  by 

V  -  f  Vv  (t)  +  (1  •  f)  VA(t)  (3.32a) 

E  =  f  Ev(t)  +  (1-f)  E  (t)  (3.32b) 

In  general  one  wishes  to  determine  P(V,E)  and  the  coupled  set  of  equations 
for  doing  this  are  Eqs.  (3.10)  (with  Eq.  (3.16)),  (3.23),  (3.28),  (3.29), 
(3.31),  (3.32a),  and  (3.32b).  These  equations  involve  the  7  variables 
t,  p(t),  v  (t),  v At)  9  f(t),  E  ( t) ,  and  E.(t). 

V  Jt  V  x 

Two  of  the  critical  parameters,  and  P^,  remain  undetermined.  In  the 
overall  equation  of  state  as  visualized,  Figure  1.2,  it  is  necessary  to 
join  the  vapor-liquid  and  solid-liquid  mixed  phase  regions  at  zero  pressure. 
Doing  this  serves  to  determine  the  critical  density  p  .  Eq.  (3.29) 
evaluated  at  the  reduced  melting  temperature  is 

p£  (V  58  pc  [L  +  c  (1  •  Si)1/3  +d  (1  *  V  ]•  (3>33) 

where  the  densities  here  are  actual  densities,  and  this  can  be  solved  for 
p  .  c  and  d  are  not  generally  known  for  the  metals  and  in  all  numerical 
work  it  has  been  assumed  that  the  values  c  ®  1.8,  d  =  1.75,  appropriate 
for  the  alkali  metals,  also  apply  for  aluminum,  titanium  and  beryllium. 

This  gives  critical  densities  of  0.557,  0.958,  0.396  gm/cc,  respectively. 
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The  beryllium  value  of  0.  J9f»  corresponds  to  Lhe  liquid  density  value  ol 
1.02  1 1  fin  *rable  2*3.  i  he  critical  pressures  are  now  tound  from 

T  -  Z  RT  /V  .  (3.34) 

c  c  c  c 

Table  1.3  summarizes  the  set  of  critical  parameters  which  have  been  deter¬ 
mined  tor  the  three  metals  considered. 


Table  3.3:  Critical  Parameters  for  Aluminum,  Titanium  and  beryllium 


Metal 

A luminum 

Titanium 

Beryl  lium 

O 

Temperature  (  K) 

8000 

1  lbOO 

8000 

Density  (gm/cc) 

0.3569 

0.9580 

0.3964 

Pressure  (kbar) 

4.  161 

5.846 

8.857 

Compressibility  Factor  (Z^)  0.307  1 

0.307 1 

0.3067 

Slope  of  V.P.  Curve 

(at)  3.775 

4.068 

3.745 

Parameter  & ^ 

-1.020 

-0.5428 

-1.  112 

Parameter  3 

5.626 

5.626 

5.638 

It  should  be  noted  that 

the  formulation  as 

described  doesn1 

1 1  give  exactly 

zero  pressure  at.  the  me 

lting  temperature. 

Values  obtained 

are  6.0  x  10  ^ 

1,9  x  10  ,  and  1.0  x  10  atm  for  aluminum,  titanium  and  beryllium, 

respectively.  Although  negligible,  these  values  have  been  subtracted 
from  the  calculated  pressures  in  the  vapor  phase  and  the  vapor- liquid 
mixed  phase  to  ensure  zero  pressure  on  the  axis.  Furthermore,  the  vapori¬ 
zation  envelope  now  terminates  at  a  very  large  but  finite  volume  rather 
than  tailing  off  to  an  infinite  volume. 
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X lie  problem  of  determining  the  critical  pressure  from  the  vapor  pressure 
equation  has  already  been  mentioned  and  by  following  the  procedure  described, 
the  *»o-called  boiling  temperatures,  instead  of  corresponding  t  o  atmos¬ 
phere,  occur  at  pressures  of  1.2,  0,32,  and  2.8  atm  for  aluminum,  titanium 
and  beryllium,  respectively.  The  boiling  temperatures  at  one  atmosphere 
are  now  2694°K  (2736),  3852°K  (3550),  and  2553°K  (2757). 

2.3  SOLID  PHASE 

An  equation  of  state  for  the  solid  phase  is  presumably  known  and  that  which 
is  required  for  the  following  formulation  is  an  equation  valid  around 
energies  where  melt  begins.  The  effect  of  any  change  in  crystalline 
structure  in  the  solid  phase  is  a  separate  problem  which  should  be  in¬ 
cluded  in  specifying  the  solid  equation  of  state.  Also,  in  the  idealized 
one -d imensiona 1  propagation  problem,  the  stress  in  a  direction  normal  to 
the  direction  of  propagation  is  different  from  that  in  the  direction  of 
propagation.  This  is  briefly  considered  in  Appendix  B.  For  the  following, 
however,  it  is  assumed  that  an  expression  for  the  hydrostatic  pressure, 
valid  at  the  higher  energies,  is  known. 

Various  forms  for  the  solid  equation  of  state  have  been  used  and,  in 
general,  they  are  of  the  form 

P(o,E)  =  cu  +  F(u.)  +  c; (Ti )  E,  (3.35) 

where  V.  =  p/p^  an^  u  38  -  -  1.  The  particular  form  assumed  for  the  purpose 
of  making  calculations  is 

P  (p ,  E)  =  Cn  +  Da2  +  Sa3  Hi  E  (3.36) 
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with  C,  l),  S  and  ti  being  constants  t  m*  a  particular  material.  Values  fur 


d  i  MIIIX  Itvllll  , 


1 —  ....  1  1  J  . 

1.  y  A  .  *  V4»l» 


A  polynomial  form  for  r(u)  ir.  unwally  adequate  for  fitting  the  data 
around  normal  densities  but  the  condition  that  the  pressure  be  a  mono- 
ton  leal  Ly  increasing  function  ot  the  density  is  often  violated.  This 
occurs  both  for  titanium  and  beryllium  where  the  constant  S  has  negative 
values.  Hence  the  P  -  a  isoenergy  curve  passes  through  a  maximum  at  a 
relatively  high  density  and  this  poses  a  problem  for  the  proposed  formu¬ 
lation.  Alternatives  are  to  alter  the  solid  equation  of  state  so  it  is 
tnonotonica 1 ly  increasing  function  of  the  density  or  to  modify  the  formu¬ 
lation.  This  problem  is  considered  again  in  the  next  section. 


As  has  already  been  mentioned,  there  is  the  condition  that  the  solid  den- 

ity  and  energy  at  which  the  material  begins  to  melt  at  zero  pressure 

must  be  compatible  with  Eq .  (3.35).  They  cannot  be  independently  specified 

and  the  suggested  procedure  is  to  specify  one  of  them  and  calculate  the 

other.  Discrepancies  introduced  by  doing  this  are  usually  not  large. 

These  quantities  will  be  designated  by  n  and  E 

1  so  so 

2.4  SOLID-LIQUID  MIXED  PHASE 

The  only  condition  which  must  be  strictly  adhered  to  in  describing  the 
solid  to  liquid  phase  transition  is  the  Clapeyron  equation  which,  in  this 
case,  is 


dP  _  AH 

JT  "  T  (V  -  V  ) 
Jl  S 


0.37) 
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£H  is  now  the  heat  of  fusion,  in  general  a  function  of  the  temperature, 
and  subscripts  s  and  i  refer  tc  the  solid  and  liquid  phases  respectively. 
In  order  to  define  the  mixed  phase  region  further,  it  is  necessary  to  make 
several  assumptions. 


It  was  pointed  out  in  Section  2.4  that,  at  the  normal  melting  temperat ure , 
the  slope  of  the  transition  curve  is  quite  large,  100  to  200  atiti/0K,  for 
the  metals  considered,  It  probably  becomes  larger  at  higher  temperatures 
but  the  assumption  is  made  that  it  remain  constant  and  equal  to  its  value 
at  zero  pressure,  Hcuiee,  the  equation  for  the  solid  to  liquid  transition 
curve  is 


P 


*  V 

o 


(’3.38) 


is  the  liquid  volume,  at  the  zero  pressure  melting  temperature,  which 
is  either  specified  or  can  be  estimated  as  some  reasonable  fraction  like 
6.57.  larger  than  the  solid  volume.  The  volume  difference  on  melting  is 
denoted  by  The  heat  of  fusion  at  zero  pressure  is  “  H  (T^) -H^ (T^) 

and  this  is  also  equal  to  the  internal  energy  difference  denoted  by  /.E  , 

A  direct  consequence  of  this  assumption  is  obtained  from  Eqs.  (3.37)  and 
(3.38),  with  H  *  E  +  PV ,  and 

E.(T)  -  E  (T)  /.E 

%  X  ...  s -  =  — - -  (  i  3  q\ 

.IV  V,  (T)  -  V  (T )  jV  *  u  ' 

i  S  O 

In  other  words,  the  ratio  of  the  energy  difference  to  the  volume  difference 
between  the  liquid  and  solid  at  the  temperature  T  is  independent  of  the 
temperature . 
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It  has  been  observed  that  ,  at  least  tor  the  lower  temperatures  near  the 

normal  melting  temperature,  both  the  melting  entropy  AH (T)/T  and  the 

16 

volume  change  .  V(l)  are  nearly  constant.  Hence,  assuming  Eq  .  ('1.3V)  is 

correct,  the  internal  energy  difference  L K(T)  is  also  nearly  constant. 


Kk<T)  =  Eso  +  C(T-Tm)  -  -  1  J  .  (3*43) 

where  1C  is  the  difference  in  specific  heats  at  T  .  Each  internal  energy 
starts  ouL  with  the  right  slope  and  becomes  linear  in  temperature  at  large 
temperatures,  In  general  AC  is  quite  small  and  hence  the  internal  energy 
difference  doesn't  change  much  for  these  assumed  forms.  Similarly  AV 
doesn't  change  much  and  the  melting  entropy  is  nearly  constant. 
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The  density  for  the  solid  on  the  melting  curve  is  obtained  by  equating 
Eqs.  (3.35)  and  (3.38)  at  the  temperature  T, 


Ps(ps(T),  Es(T))  =  P(T). 

0.44} 

In  the  melting  region  the  volume  V  and  internal  energy  f 

through  the  liquid  fraction  f  by 

are  related 

V  =  f  v£  +  <i-Ova, 

(3.45a) 

E  >  f  +  (1-f)  Eb. 

(3.45b) 

All  quantities  on  the  right  of  these  equations  are  functions  of  the  tempera¬ 
ture  T. 

The  coupled  set  of  equations  for  finding  P (V  f  E)  in  this  mixed  phase  region 
ave  now  Eqs.  (3.38),  (3.39),  (3.42),  (3.43),  (3.44),  (3.43a),  and  (3.45b), 
involving  the  7  variables  T,  P(T),  V.C'i'),  V  (T),  t  (j'),  E,(T),  and  E,  (T)  . 

S  jj,  S  Jt 

Again,  the  simplest  iterative  procedure  converges  rather  rapidly. 

Qualitatively,  this  formulation  has  the  effect  of  bending  the  transition 
region  toward  higher  densities  at  the  higher  pressures.  One  result  of  this 
is  that  a  phase  transition  can  occur  at  normal  density  for  energies  some¬ 
what  higher  than  the  zero  pressure  values  and  rather  high  pressures,  This 
is  illustrated  in  Table  3.4. 
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i  :l  I)  I  i  '  \n  I  himij  i  t  v  M*i  I  4  n  r»  l-’nijr  [ii  p  o  fjnH  Pr^oaiM'pe 


Me  t  a  i 

A  i  urn  nun 

11  tan J urn 

He ry 1 11 urn 

t.iu-rgy  at  Zero  Pressure  Kal'gm) 

11  1 3 

280 

878 

Ijiergv  at  Normal  Density  (cal/gm) 

340 

42  7 

1360 

Pressure  at  Normal  Density  (kbar) 

81 

87 

120 

temperature  at  Normal  Density  (  K) 

1440 

2820 

2170 

In  Section  3.2  the  reference  energy  for  the  vapor  was  undetermined  and  it 

is  now  chosen  so  that  Lhe  liquid  energy  is  equal  to  1C,  ,  all  energies 

k  o 

being  referenced  to  zero  energy  for  the  solid  at  normal  density  and  zero 
pressure . 

hushing  this  model  still  iurther,  to  higher  densities,  pressures,  and 
energies,  requires  that  the  solid  equation  of  state  be  assumed  correct 
under  these  extreme  conditions.  In  doing  this  the  problem  mentioned  in 
the  last  section  arises  and  if  the  solid  isoenergy  curve  is  not  a  mono¬ 
tonies  lly  increasing  function  of  the  density  one  eventually  reaches  a 
situation  where  Llq .  (3.43)  no  longer  has  a  solution.  This  limitation  is 
probably  not  serious  for  the  solid  phase  end  the  solid- liquid  mixed  phase, 
where  the  densities  are  not  expected  to  be  too  large,  but  it  does  prevent 
one  from  determining  the  boundary  of  the  liquid  phase  for  high  energies. 

Aside  from  this  problem,  the  proposed  model  completely  describes  the 
solid-liquid  transition  region  and,  in  particular,  fixes  the  pressure, 
density  and  energy  on  the  liquid  side  of  this  region. 


3-18 


2.b  LIQUID  PHASE 


As  has  been  mentioned,  the  liquid  phase  seems  to  be  the  most  difficult 
phase  to  describe*  In  principle,  the  equation  of  state  formulation 
developed  to  this  point  defines  the  boundaries  of  the  liquid  phase 
region  and  determines  the  thermodynamic  variables  on  these  boundaries. 
Any  assumption  made  for  bridging  the  gap  is  highly  questionable,  however 

A  possible  approach  is  to  construct  isotherms  through  the  liquid  region 
and  find  the  energies  along  the  isotherms  from  Eq .  (3.8)  .  The  proposed 
formulation,  however,  determines  the  energies  at  the  end  points,  along 
the  liquid  sides  of  the  mixed  phase  regions,  and  obtaining  agreement  by 
this  approach  may  be  rather  formidable.  Since  the  quantity  of  primary 
interest  is  internal  energy,  a  simpler  procedure  is  to  develop  empirical 
expressions  for  curves  along  which  it  is  constant. 

Various  expressions  have  been  considered  and,  since  the  choice  is  some¬ 
what  arbitrary,  there  is  a  tendency  toward  simplicity  as  long  as  the 
results  are  "reasonable." 


One  simple  expression  assumes  that  the  pressure  is  a  linear  function  of 


the  density  for  a  given  energy.  Hence,  denoting  this  pressure  by  P^(p,E), 


P,(P,E)  =  — ^ ^  p 
pJlM  P£B 


3  P 
xH  IB 


P£M  °£B 


0.46) 


where  the  subscripts  M  and  B  denote  the  melting  and  boiling  sides  of  the 
liquid  region,  respectively.  All  subscripted  variables  are  functions  of 
the  energy  E.  For  energies  greater  than  that  at  the  critical  point, 
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.1  *  j  and  V  >  P  .  Pressures  obtained  with  this  equation  seem  quite 

<13  c  ili  c 

high,  particularly  at  the  higher 

A  second  simple  expression  results  from  taking  the  logaiillmi  of  the 
pressure  as  a  linear  function  of  the  logarithm  of  the  density.  Denoting 


this  pressure  by  P9(o»K)> 

in  (  W 

«n  P>’E)  =  7n  (TV)  *"  0  - 

V£M 


* 11  °jB  in  P£M  '  fn  g£M  iXl  P£B 


£n  (o£M/p£B) 


(3.47) 


This  gives  rather  reasonable  pressures  for  the  higher  energies  but  leads 
to  a  crossing  of  curves  of  constant  energy  for  the  lower  energies.  In 
fact,  the  curve  may  cut  down  through  the  melting  region. 

A  third  expression,  and  the  one  included  in  the  equation  of  state  program, 
is  obtained  by  taking  a  linear  combination  of  the  logarithms  of  the 
presures  P^  and  .  Hence 

in  P  (p,E)  =  F (E)  £n  P^E)  +  (l-F(E))  in  P2  (p,E)  (3.48) 

where  F(E)  is  somewhat  arbitrary.  Reasonable  results  have  been  obtained 
by  limiting  it  to  the  range  0  to  1  and  determining  it  from  the  condition 


iP(p.E) 


JrcJfL 

.(E)  -  P., 


(3.49) 


By  choosing  it  this  way,  the  curves  of  constant  energy  do  not  cross  and, 
in  general,  Eq .  (3.46)  is  predominant  at  low  energies  while  Eq .  (3.46) 
determines  the  pressure  at  high  energies. 
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SECTION  4 

EQUATION  OF  STATE  SUBROUTINE 


An  equation  of  state  subroutine  for  finding  the  pressure  P  for  a  given 
volume  V  and  internal  energy  E,  according  to  the  formulation  described 
in  the  previous  section,  has  been  programmed  and  checked  by  using  it  in 
two,  rather  simple,  main  programs*  It  has  not  been  tested  in  a  PUFF 
hydrodynamic  program.  In  writing  this  program  the  emphasis  was  placed  on 
obtaining  a  workable  program  rather  than  on  efficiency  and  hence,  there 
is  undoubtedly  room  for  improvement.  A  listing  of  the  subroutine  is 
given  in  Appendix  C. 

The  required  input  information  is  given  in  Table  4.1.  The  equations 
which  are  solved  are  scattered  throughout  Section  3  with  jus t if icat ion 
for  the  assumed  forms.  They  are  assembled  here  for  clarity.  The 
density  p  and  volume  V  are  both  used,  depending  on  which  is  more  con¬ 
venient,  and  one  is  just  the  reciprocal  of  the  other.  In  some  instances 
the  reduced  values  are  also  used.  A  consistent  set  of  units  is  implied. 
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Table  4.1:  Input  In formation  for  Equation  of  State  Subroutine. 

Notation 


Description 

L’nits  Formulation 

Subrou  tii 

Atomic  Weight 

gm/mo  le 

M 

WT 

Normal  Density 

gm/cc 

DO 

Critical  Temperature 

°K 

T 

c 

TCK 

Coefficient  in  Liquid  Density 

Eqnat ion 

None 

c 

Cl 

Coefficient  in  Liquid  Density 

Equat ion 

None 

d 

Dl 

(  fl  } 

Parameter  in  Vapor  EOS 

None 

k 

0 

KO 

Melting  Temperature 

°K 

lM 

TMK 

Solid  Enthalpy  at 

kca 1/mole 

",<V 

HSM 

Liquid  Enthalpy  at 

kcal/mole 

w 

nui 

Vapor  Enthalpy  at  T^ 

kca 1/mo le 

W 

HVN 

(b) 

Solid  Density  at  T^ 

gm/cc 

P/V 

DSM 

( c ) 

Liquid  Density  at 

gm/cc 

p£<TH^ 

DLM 

Boiling  Temperature 

°K 

TBK 

Liquid  Enthalpy  at 

kcal /mo le 

W 

HLB 

Vapor  Enthalpy  at  T 

D 

kca  1/mole 

W 

HVB 

(d) 

Coefficient  in  Solid  EOS 

kbar 

c 

C 

(d) 

Coefficient  in  Solid  EOS 

kbar 

D 

D 

Coefficient  in  Solid  EOS 

kbar 

S 

S 

Coefficient  in  Solid  EOS ^  ^ 

kbar/ca  1/gm 

G 

c 

Zero  (Make  initial  ca leu lat ions) 

None 

INI 

Legend  for  Table  4*1. 

a.  If  this  is  zero,  it  is  set  equal  to  3. 

b.  If  this  is  zero,  it  is  determined  from  the  solid  EOS  with  E  '  H  (T.J 

s  M 

If  a  value  is  specified  the  corresponding  energy  is  determined 
from  the  solid  EOS  and  .ill  other  energies  are  renormalized. 

c.  If  this  is  zero,  it  is  taken  to  be  6.57,  less  than  the  solid 
density. 

d.  The  solid  EOS  is  assumed  to  be  of  the  form  given  in  Eq.  (3*36). 


Figuve  4.1  shows  the  various  regions  in  a  P-V  diagram  and  defines 
certain  quantities  useful  in  specifying  the  boundaries  of  these  regions. 
The  complete  set  of  equations  is  as  follows: 

^r^surejBCalcu^aXion£ 

Solid  Phase:  71  *  p /p  ,  u  =  Tl  -  1. 

Pg  (p,E)  =  Cu  +  Du2  +  Su3  +  GT|E.  (4.1) 

Determine  P  =  Ps  (p,E) 
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Se  lid -Liquid  Mis?**  ■ 


Pm  •  aR  /aV  (T/T. .  -  1) 

”  o  “  o  M 


(4.2a) 


Pg(Ps(T),  ES(T))  -  P (T ) 


(4.2b) 


E„(T)  -  E  (T)  /£ 

- a - o_ 


V  (T)  -  V  (T) 

l  S  O 


AC  T, 


Ei(T)  "  E£o  +  C  VT/Vl)  *  2 


(Tm/T-1) 


(4.2c) 


(4. 2d) 


AC  Tm 

Es(T)  =  Eso  +  C  1m('1/tm-1>  +  — 7  VT  *  1} 


v  =  fe(T)  VT)  + 


l-f£(T)jVs(T) 


(4.2e) 

(4.2f) 


E  *  f£(T)  E£(T)  +[l-£t  (T )jEg(T) 


(4.2») 


Determine  T,  V } ,  V*,  E  ,  Eg,  f£,  P  -  P(T). 


Liquid  Phase: 


Pl(p,E)  =  -V~  V  p  -  to.  1  p^i  F^b 


j2M 

A(E)p  -  B(E) 


piM  '  P£B 


„  ,  pv _ ~IM'  ~IB  ’  .  P^B*  IM 

i,n  P9  (o  > E)  -  _  _  j  ,  x  Jtn  p  •  n\  t  7  T 

-  £  p  £M/p  £B  £  (pi>l//pXB) 


(4.3a) 


“  PoB£n  Pi>M  “  ■£n  P„M  in  PfB 


(4. 


=  C(E)  in  p-D(E) 


3b) 
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F(F>  -T-iH  C(E)  -  -  i/r  -  C(E)  -  A(E)  j  (A. 3c) 

•'.'M  '  M  C  lo  J'  •'  )M 

M  P,  (p,E)  =  l-(i:)  in  P.(o,E)  +f  l-F(li)  i  £n  P„(C,E)  (A. Ad) 

l  1  _  J  ^ 

Determine  P  =  P,(o , K) 

Liquid-Vapor  Mixed  Phase:  (Reduced  densities) 

p(c)  =  exp|  (or-u  |)  ( 1  -  1/t )  +  Or |  Cit  cj  (4.4a) 

p  (o  (t),t)  ®  p(t)  determines  a  (t)  (4.4b) 

.  .u)  =  1  +  c  ( 1  "t- )  +  d(l-t)  (A. Ac) 

F.  (t)  =  12  +  C  T  t  -  Z  RT  T  (k  +  2k. /t)  p  (t)  -  k.  p2(t)/ti  (A. Ad) 

v  o  vc  cc  o  1  Mv  2  v  J 

Ee(t)  =  Ev(t)  -  ZcR’Lc  rvyU)  -  vi(t)j  p  (t)  k-ttl)/t  +ax  -  lj  (A.Ae) 
V  =  £v(C)  Vv(t)  +  f  l-fv(t)jV£(t)  (A  .  C) 

E  -  fv(t>  Ey(t)  +[l-£v  (t)j  Ee(t)  (A.4g) 

Determine  t,  vy,  v£,  Ev>  E^,  fy  P(V,t)  =  Pc  [p(t)  -  pvoj 
Vapor  -Phase  :  (Reduced  density') 

P  (?,t)  *  p.t/Z  (1-bp  +  b '  ^ 2 )  -  (k  +  k./t)  r2  -  k  (t-l/t)p3  (A. 5a) 

v  c  o  1  L 


E  -  E0  +  CvTct  -  Zc  RTc  [<kQ  +  2^/t)  p  -  k2  p2/tj 

Determine  t,  P  a  P  |  p  (D,t)  -  p  i. 

C  L  v  rvoJ 

Calculations  of  Boundary  Variables 
Solid-Liquid  Mixed  Phase  Region: 


Solve  Eqs,  (4.2a),  (4.2b)  and  (4.2e)  with  Eg  =  E,  for 

VE)  =  Vs(T)‘ 


E  <  E  v- 
so 


E 


i  o 


V 

max 


(E 


E  )  V „ 
so  %  o 


(E  -E)  V 
lo  so 


(E 


io 


-  E  ) 
so 


If  V  >  V  (E)  ,  P  5  0. 
max 


E.  <  E 
to 


Solve  Eqs.  (4 . 2a)  - (4 . 2e )  with  *  E,  for  V^(E)  - 

VH  -  KT). 


Liquid -Vapor  Mixed  Phase  Region: 


E.  <  E  <  E 
to  c 


Solve  Eqs.  (4 .4a) -(4 .4e)  with  E  =  E,  for  V^(E)  * 


p,B(E)  -  Pc  !  p(t)  -  Pvo 


.5b) 


(I), 


/o  ,  (  t  )  , 
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/  >  y  (r  -  e)v  !  -p*  ) 

max  "  t'o  vo  Vo  *  *’ oj  /  '“vo  ‘Ho 


It'  V  '»  V  ( I’)  ,  r  =  0. 

in  ax 


K  v  K 
e 


V ^ |  < K )  ~  V  •  Solve  Eqs.  (4.5a)  and  (4.5b),  with  p  =  1, 


tor  VK)  =  PcLPv(I,t)  ‘  PvoJ 


V ,  vT  v  V 

;o  c 


Solve  Eqs.  (4 . 4a)  -  (4 , 4o)  ,  with  p  =  V  /V,  for  I£  (V)  ~  K  .(C) 

Jr.  C  {.  13  / 


V  <  V  •:  V 

C  VO 


Solve  Eqs.  (4.4a),  (4.4b),  and  (4.4d),  with  n  =  V  /V,  for 

w  v  c 

iyv)  =  Ev(t). 


Calculations  of  Equation  of  State  Parameters. 


If  o  0\«)  <  0;  E  =  H  (T  )  ,  P  (p  ,  E  )  =  0  determines  p 

ws  M  -  so  s  M  s  *so  so  Kso 

v  0;  p  *  p  (1‘  )  ,  P  (p  ,  E  )  =  0  determines  E 

’  so  rs  M  sHso  so  so 


E ,  *  E  +  11,(1*  )  -  II,  Cl*) 

to  so  jC  M  s  M 


If  p  ('!’  )  o  0;  =  0.935  p  .  (A  different  fraction  could  be  a  sod.) 

<  M  -  ;o  s  o 

*  0;  =  o,dM) 

a.o  M 


+  “  “Vj 
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*HB  “  W  *  W 


<*R  +■  0.0838  (OrR  -  3.75)  (36/tg  +  6tR  -  42)  *  AHfi/KTc  deter¬ 


mines  a 


zc  =  \i.n  +  0.26  (*R  -  7)  | 


-1 


P  -  Z  KT  3  /M 
r  e  c  r  c 


cp  ~  [UV(1B')  ‘  hv-1m^J  <VV*  cv  “  cp  "  R 

ctM[W  'Wj  <VV 


AC  -  Cp  ■  ^ 


Of  *  +  (Tc-TR)j  RTC,  orL  =  #.C/R 

2  *  5(3  3-l)/(:+l)J  determines  5 


b  *  (3  3  -  60  -  l)/3(3s-l),  b 1  -  (--3)/(3gt-l) 


If  k  =  0,  k  ■  B*  If  k  i  0,  k  =  k  . 
o  o  o  o  o 


kl  °  fe  "  lV  k2  *  2  (1  *  ko  +  23  ’  0f) 


Let  E  =°  H  (T  )  -  CL  Tt,  =  E 1 
o  v  B  P  B  o 


Solve  Eqs .  (4 . 3a) - (4. 3e) >  with  t  =  t*4,  for  p  -  ?  o  (t),  p 

M  vo  c  rvo 

E  =  E'  +  [~E  -  E,(t)  ,  E  =  (E  -E  )  +  E  (t). 

o  o  [  fLo  *  j  vo  o  o  v 


p(t ) 
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Hi  is  re  m»nna  l  i  &ai  i  on  of  t'Ut'itv  assures  that  K  “  K  (t  ) 

A*U  i  N 

E  s  E  +  Ct  'i  -  2  R  T  Ik  +  2k.  -  k-  i 
c  o  V  c  c  c  l  o  1  2  J 


•:  *  e.  -  k 

o  so 


*  V  -  V .  ‘V 
O  «.!)  so 


C .  55  H  (  i*  )/('!\.  -  2  98) 
s  s  M  N 


In  essence  the  subroutine  does  three  things: 

1. )  The  first  Line  through  it  calculates  the  critical  para¬ 
meters,  and  other  parameters  appearing  in  the  equations,  from 
a  limited  amount  of  input  information.  (This  is  only  done 
once  and  these  calculations  could  be  placed  in  the  main  pro¬ 
gram.) 

2. )  It  calculates  the  necessary  boundaries  of  the  various 
regions  of  the  phase  diagram  for  the  given  values  of  the 
energy  E  and  volume  V. 

3).  It  determines  which  phase  or  mixed-phase  region  the 
values  of  the  volume  and  energy  fall  in  and  calculates  the 
pressure  accordingly. 

The  checking  starts  with  the  solid  phase  and  works  toward  the  vapor 
phase.  Hence  a  minimum  amount  of  time  is  wasted  in  determining  a  pressure 
in  the  solid  phase. 
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As  written,  the  subroutine  give  the  pressure  in  kbar  tor  the  volume  in 
cc/gm  and  the  energy  in  e « 1  / gtu .  Thio  la  obtained  by  the  statement 

•Ties*, ui’fc"  -  r(”Vn lumc",  "Energy")  (4 . 1  > 

where  any  designat i on  for  "Pressure",  "Volume",  and  "Knergy"  in  the  main 
program  may  be  used. 

The  explicit  form  for  the  solid  equation  of  state  given  by  Hq .  (4.1)  has 

been  considered  but  the  more  general  form  of  hq .  0.35)  only  requires 
reprogramming  the  functions  F(fi  ),  ti(r),  and  their  derivatives.  These 
are  designated  by  IT  (MU)  ,  GG(ETA),  FFP(MU) ,  and  GCP(ETA)  respectively. 
Corresponding  changes  must  also  be  made  in  the  input. 
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SECTION  5 

NUMERICAL  RESULTS 


The  equation  of  state  subroutine  described  in  Section  4  and  lisln!  in 

Appendix  C  has  been  used  in  two  programs  for  the  purposes  of  checking  it 

out  and  obtaining  some  numerical  results  for  aluminum,  titanium,  and 

beryllium.  In  the  process  of  doing  this,  certain  parameters  appearing  in 

the  formulation  were  extracted.  These  are  given  in  Table  5.1.  The  critical 

constants  for  aluminum  from  Ref.  2  are  T  =  6000°K,  p  *  0*590  gm/cc, 

c  c 

P  3  3.11  kbar,  and  Z  *  0.292.  In  view  of  the  lack  of  data,  the  differences 
c  c 

are  probably  not  significant. 

Figures  5.1  -  5.3  show  curves  of  constant  energy  for  aluminum,  titanium, 
and  beryllium,  respectively.  The  -n  P  vs  ^n  P  plots  best  show  the  details 
of  the  liquid-vapor  mixed  phase  and  the  vapor  phase  where  the  pressures  and 
densities  are  relatively  low. 

The  high  pressure,  high  density  region  for  aluminum  is  shown  in  Figure  5.4. 

On  this  plot  the  liquid-vapor  mixed  phase  is  completely  lost.  In  addition 
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iau  m  j  •  x  . 


odicu  laieu  rj^uci  u  .lou  ui  olulc  rdiaumut  s 


Metal 

Aluminum 

Titanium 

Beryllium 

Critical  Parameters: 

Temperature  (°K) 

8000 

11600 

8000 

Density  (gm/cc) 

0.5569 

0.9580 

0.3964 

Pressure  (kbar) 

4.161 

5.846 

8.857 

Compressibility  Factor  (Z^) 

0.3071 

0.3071 

0.3067 

Internal  Energy  (cal/gm) 

3016 

2802 

8899 

Slope  of  V.P.  Curve  (pt) 

3.775 

4.068 

3.745 

Parameter  3 

5.626 

5.626 

5.638 

Zero  Pressure  Parameters: 

Solid  Density  (gm/cc) 

2.537 

4.251 

1.732 

Liquid  Density  (gm/cc) 

2.380 

3.974 

1.619 

Vapor  Density  (gm/cc) 

2.124  x  10' 15 

5.682  x  10' 10 

7.079  x  10' 

Solid  Energy  (cal/gm) 

212.6 

279.5 

878.3 

Liquid  Energy  (cal/gm) 

307.0 

356.9 

1189 

Vapor  Energy  (cal/gm) 

2994 

2451 

9032 
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PRESSURE  -  KB 


FIGURE  5.1  ALUMINUM  ISOENERGY  LINES.  PARAMETER  IS  ENERGY 
IN  CAL/GM. 
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PRESSURE 


FIGURE  5.3  BERYLLIUM  ISOENERGY  LINES.  PARAMETER  IS 
ENERGY  IN  CAL/GM. 
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PRESSURE  (Mb) 


/  .no 


FTCURK  r>.4  ALUMINUM  ISOENERGY  LINES  AND  AD  I A  BATS 


5- 


to  curves  of  constant  energy,  calculated  adiabats  are  also  shown.  At  the 
nigh  pressures,  considerable  PdV  work  is  done  and  hence  there  is  a  marked 
change  in  internal  energy  along  an  adiabat. 

Figure  5.5  shows  adiabat9  for  aluminum  which  cut  down  through  the  liquid- 
vapor  mixed  phase  region.  At  these  low  pressures  there  is  little  distinction 
between  adiabats  and  curves  of  constant  energy.  Provided  the  formulation 
is  reasonably  correct,  there  arc  marked  breaks  in  the  curves  when  vaporiza¬ 
tion  begins . 

Figure  5.6  is  another  plot  of  the  adiabats  of  Figure  5.4  where  the  internal 
energy  is  shown  as  a  function  of  density.  It  more  cl.early  indicates  the 
effect  of  vaporization  in  an  adiabatic  process.  Although  the  calculations 
at  very  low  densities  were  not  made,  It.  is  expected  that  energy  recovery 
proceeds  down  to  the  liquid  internal  energy. 

It  should  be  noted  that  the  pressures  of  Figure  5.4,  for  a  given  energy, 
are  about  a  factor  of  two  higher  than  those  obtained  in  Ref.  2.  This  is 
illustrated  in  Figure  5.7  which  gives  the  pressure  at  normal  density  as  a 
function  of  internal  energy.  Also  shown  is  the  normal  density  pressure 
obtained  from  the  solid  equation  of  state.  This  is  the  pressure  that  most 
older  versions  of  PUFF  would  predict  for  normal  density  aluminum.  The  dif¬ 
ference  is  due  to  the  empirical  treatments  of  the  so-called  liquid  phase 
and,  at  present,  there  seems  to  be  no  data  to  support  or  refute  either 
formulation . 

The  equation  of  state  program  as  written  is  not  an  efficient  program  for 
determining  curves  of  constant  energy  since  it  starts  from  ’’scratch" 
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ALUMINUM  AD IA BATS  SHOWING  I  RGY  RECOVERY 


for  each  pair  of  values  of  the  volume  and  energy.  If  one  wished  to  repre¬ 
sent  the  of  ststs  in  tibuls?  forts,  ss  in  ®*f f  2,  2  scii f isi  pro¬ 

gram  for  obtaining  these  tables  should  be  prepared* 

The  adiabats  were  calculated  according  to  the  difference  equation  connect¬ 
ing  points  1  and  2 , 

P  +  p 

E2  -  Ex  +  —2 - -  (V2  -  vp  -  0  (5.1) 

where  the  *  native  process  was  written  into  the  main  program. 

All  calculations  were  made  on  the  Aeronutronic  Time-Sharing  System  which 
utilizes  the  GE  420  as  the  main-site  computer.  A  rough  estimate  of  the 

average  time  involved  to  obtain  one  pressure  for  a  given  volume  and  energy 

-4 

is  0.25  sec,  The  accuracy  imposed  on  iterations  in  the  subroutine  was  10  1, 
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SECTION  6 

SUMMARY 


In  summary,  two  main  objectives  of  the  present  investigation  have  been 
achieved : 

1)  An  analytical  model  of  an  equation  of  state,  with  phase  changes, 
for  pure  metals  has  been  developed. 

2)  An  equation  of  state  subroutine  which  yields  the  pressure  for  a 
given  specific  volume  and  energy  has  been  written. 

In  developing  the  equation  of  state  model,  the  effects  of  dissociation  and 
the  possibility  of  sublimation  have  been  neglected.  Hence,  at  best,  it  can 
only  apply  for  pure  materials  which  normally  melt  and  whose  vapor  consists 
of  atomic  species.  The  lack  of  data  for  high  temperatures  and  pressures 
necessitated  many  assumptions  in  order  to  complete  the  model  and  in  most 
cases,  plausibility  is  their  only  justification. 

The  data  required  to  determine  certain  parameters  appearing  in  the  formula¬ 
tion  are  essentially  no  different  from  those  which  have  been  used  in  the 
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past.  Now,  however,  they  are  used  to  describe,  in  a  reasonable  fashion, 
the  phase  changes  from  solid  to  liquid  and  liquid  to  vapor. 

The  only  metala  considered  to  date  have  been  aluminum,  titanium  and  beryl¬ 
lium.  While  no  difficulties  arose  for  these  metals  there  Is  no  assurance 
that  there  isn't  some  material  which  will  tax  the  proposed  model.  A  few 
preliminary  calculations  for  a  new  material  are  helpful  to  check  for 
"reasonability." 

Both  the  equation  of  state  model  and  the  computer  subroutine  evolved  over 
a  period  of  time.  As  a  result,  the  program  isn't  as  "clean"  as  it  could 
be.  The  main  effort  was  to  obtain  a  workable  subroutine  and,  hence,  many 
of  the  iterative  procedures  can  undoubtedly  be  improved. 

The  numerical  results  of  Section  5  illustrate  the  type  of  information  that 
one  can  obtain  with  this  program.  It  is  noted,  however,  that  the  program, 
as  written,  gives  the  pressure  for  a  given  volume  and  energy  with  no  use 
being  made  of  any  prior  information  concerning  the  location  of  the  previous 
point  in  the  phase  plane.  Therefore,  it  is  quite  inefficient,  e.g.,  for 
obtaining  curves  of  constant  energy. 
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APPENDIX  A 


CRITICAL  POINT  MEASUREMENTS  USING  INERTIAL  CONFINEMENT 


Above  about  2000°C  no  equation  of  state  measurements  in  the  liquid  vapor 
region  have  been  made  using  conventional  techniques  because  of  the  lack  of 
suitable  refractory  containers  and  instrumentation.  The  critical  tempera¬ 
ture  of  many  metals  is  well  above  this  conventional  technique  limit.  How¬ 
ever,  it  appears  feasible  to  make  measurements  at  these  higher  temperatures 
and  pressures  by  use  of  inertial  confinement  for  times  of  the  order  of 
microseconds.  A  number  of  requirements  must  be  met  if  this  approach  is  to 
succeed,  the  principal  one  being  a  means  of  uniformly  heating  the  test 
material  in  a  time  short  compared  to  the  inertial  confinement  time.  There 
\e  a  number  of  ways  in  which  this  might  be  done. 

With  such  a  means  of  rapidly  heating  a  confined  material  of  less  than  normal 
density,  the  material  may  be  driven  almost  instantaneously  to  a  pressure  above 
the  critical  pressure.  The  material  will  then  expand,  compressing  the  walls 
containing  it.  This  expansion  will  be  adiabatic  if  energy  losses  can  be 


A-l 


neglected,  so  the  pressure-density  history  will  follow  an  adiabat,  similar 
to  one  of  those  shown  for  aluminum  in  Figure  5.6.  The  pressure  history  may 
be  observed,  and  from  it  the  density  history  determined.  The  discontinuity 
Ins  Tope  of  tHe  s  e  observed  adiabats  at  the  liquid- vapor  saturation  envelope 
serve  to  locate  the  envelope,  and  from  it  the  critical  point  can  be  found. 

Specifically,  a  method  by  which  the  liquid-vapor  region  equation  of  state 
data  may  be  obtained,  utilizing  inertial  confinement,  involves  the  follow¬ 
ing  steps: 

(1)  Confine  a  thin  sheet  of  the  material  to  be  investigated  between 
two  solid  plates  of  material  of  sufficient  size  to  make  the 
geometry  effectively  one  dimensional. 

(2)  Very  rapidly  heat  the  confined  material  uniformly  throughout  its 
volume  to  generate  the  required  pressure. 

(3)  Minimize  container  wall  heating  to  reduce  background  pressures, 
and  to  allow  use  of  known  cold  properties  of  wall  material. 

(4)  Measure  or  otherwise  determine  the  specific  thermal  energy  de¬ 
livered  initially  to  the  test  material. 

(5)  Observe  the  stress  wave  transmitted  to  the  container  wall  utiliz¬ 
ing  a  suitable  fast  response  pressure  transducer. 

(6)  Deduce  the  pressure  and  density  history  of  the  test  material 
using  the  known  dynamic  equation  of  state  of  the  wall  material. 

(7)  Determine  the  thermal  history  of  the  test  material  from  the  known 
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initial  specific  energy  and  the  approximately  adiabatic  expansion 
(corrected  for  small  energy  losses). 

(6)  Locate  Che  liquid-vapor  saturation  envelope  as  breaks  in  the 
pressure-density  traces  (adiabats)  observed. 

(9)  Determine  the  critical  point  by  varying  initial  conditions  to 
trace  out  the  complete  saturation  envelope. 

Of  particular  intwi-w^^  the  present  time  is  critir;al-rcgion  data  on  metals, 
which  have  estimated  critical  pressures  up  to  about  10  Kb,  and  critical 
temperatures  up  to  about  10,000°C  (e.g.,  Figures  5.1  to  5.3).  The  critical 
densities  are  less  than  half  of  the  solid  density.  To  attain  these  pressure 
and  temperature  levels,  the  almost  instantaneous  specific  energy  deposited 
throughout  the  volume  of  material  must  reach  levels  of  3000  cal/gm  or  more. 
Two  methods  of  heating  materials  this  much,  this  fast  are  certain  to  work. 

The  first  is  by  the  use  of  radiation  from  a  nuclear  explosion,  and  the 
second  is  by  the  use  of  a  Q-switched  laser  pulse,  suitably  focused  down. 

The  specific  energy  levels  may  also  be  attained  by  dumping  a  capacitor 
(pulsed  power)  bank  and  heating  the  material  electrically  by  use  of  electron 
beam  pulses,  and  by  use  of  very  high  velocity,  very  thin  flyer  plates. 

Container  walls  must  be  of  the  order  of  a  few  millimeters  thick  to  guarantee 
inertial  confinement  (i.e.,  no  reflected  relief  stress  wave  from  outside) 
for  one  microsecond.  In  one-dimensional  geometry,  the  minimum  diameter  of 
the  heated  area  must  be  of  the  order  of  a  centimeter  to  avoid  edge  effects 
at  the  center  for  one  microsecond.  Also  required  is  a  pressure  transducer 


which  can  measure  many  ki lobar  pressures  over  a  small  point  only.  The  laser 
interferometer  transducer  recently  developed  at  Sandia^  has  this  capability. 
Others,  such  as  manganan  wire  transducers  may  also  be  adaptable  to  meet  this 
requirement . 

The  specific  energy  deposited  may  either  be  calculated  from  a  knowledge  of 
the  amount  of  energy  input  corrected  for  losses,  or  be  measured  directly 
by  some  form  of  calorimetry  or  spectroscopy.  It  is  important  that  the  energy 
L»u  uniformly  introduced  through  the  volume  of  the  material,  as  well  as  ac¬ 
curately  known  or  measured.  For  penetrating  nuclear  radiation  no  great  pro¬ 
blem  is  expected,  but  for  lasers  this  requirement  may  prove  difficult.  Un¬ 
iformity  of  energy  deposition  in  either  case  may  be  enhanced  by  using  pow¬ 
dered  target  material  to  utilize  statistical  shine- through  and  multiple 
scattering . 

To  make  quantitative  estimates  of  feasibility,  a  particular  experiment  was 
considered  using  lasers  to  heat  aluminum.  The  pressure  transducer  chosen 
was  a  laser  interferometer,  and  the  energy  delivered  was  determined  from  a 
knowledge  of  the  input  pulse  magnitude  and  energy  retention  efficiency. 

The  test  material  mass  is  limited  by  the  amount  of  energy  available  in  a 
single  pulse  of  a  pulsed  laser.  A  typical  Q-switched  ruby  laser  can  gene¬ 
rate  a  pulse  about  50  nanoseconds  long,  having  an  energy  of  about  1  joule, 
or  1/4  calorie.  This  limits  the  mass  of  material  which  may  be  heated  to 
3000  cal/gm  levels  to  about  0.1  milligram.  For  a  given  desired  initial 
density,  this  in  turn  limits  the  test  material  volume. 

The  dimensions  of  the  test  volume  is  determined  by  the  interplay  of  three 
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factors.  The  diameter  limits  the  confinement  time  before  edge  effects  be- 
gin.  The  initio  1  thickness  determine*  £li*s  mle  of  pressure  drop  due  to 
compression  of  the  chamber  wells*  Together*  the  diameter  and  thickness 
are  pinned  to  the  test  volume  limit. 

Considering  these  factors  simultaneously  gave  the  following  values  for  a 
1  gm/cc  aluminum  target  material  between  fused  quartz  chamber  walls: 

-4 

Thickness:  9  x  10  cm  0.4  mils) 

Diameter:  0.4  cm 

Confinement  time:  300  nanoseconds 

To  fix  these  values,  it  was  required  that  in  the  confinement  time  the 
volume  of  the  test  chamber  no  more  than  double. 

For  these  conditions,  the  pressure  pulse  profile  generated  in  the  fused 
quartz  chamber  walls  is  easily  within  the  measuring  capability  of  a  laser 
velocity  transducer.  The  radiation  and  conduction  losses  during  the  con¬ 
finement  time  are  of  the  order  of  5  percent,  and  thus  the  expansion  is 
nearly  adiabatic. 

The  major  problem  which  must  be  solved  if  a  laser  pulse  is  used  to  do  the 

heating  is  to  assure  uniform  heating  of  the  test  material.  This  appears 

to  be  possible  in  the  direction  of  the  laser  beam  by  using  powdered  aluminum 

particles  of  appropriate  size  (of  the  order  of  microns).  Based  on  turbid 
18 

medium  theory,  target  absorptions  approaching  50  percent  are  possible  even 
for  highly  reflecting  particles  because  of  multiple  scattering,  and  at  the 
same  time  the  distribution  in  depth  is  kept  reasonably  uniform. 
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The  use  of  a  pulsed  laser  to  obtain  valuable  equation-of-state  data  near 
the  critical  poirt  appears  feasible.  This  technique  may  be  traded-off 
against  other  approaches  to  find  the  most  promising  one*  Obtaining  equa¬ 
tion-of-state  data  by  one  of  the  inertial  confinement  approaches  will  make 
possible  location  of  critical  points  for  metals  and  refractory  materials, 
and  clarify  the  form  of  their  equations  of  state  in  the  liquid  and  liquid- 
vapor  regions. 
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APPENDIX  B 

RELATION  BETWEEN  STRESS  AND  PRESSURE 


In  general,  the  value  of  Poisson’s  ratio  v  for  a  solid  is  different  from 
1/2  and  hence,  in  an  idealized  one-dimensional  propagation  problem,  the 
stresses  normal  to  and  parallel  to  the  direction  of  propagation  are  dif¬ 
ferent*  This  is  not  the  case  for  liquids  and  vapors  where  v  =  1/2.  Since 
the  material  motion  is  governed  by  the  gradient  of  the  stress  in  the  direction 
of  propagation,  it  is  desirable  to  obtain  a  relation  between  this  stress  and 
the  so-called  hydrostatic  pressure  used  throughout  the  equation  of  state 
formulation. 


The  theory  of  thermal  stress  analysis,  in  the  limit  of  small  strains,  is 
well  known  and,  for  a  strain  s  only  in  one  direction,  the  stresses 

parallel  (c)  and  normal  (p  )  to  this  direction  are 


a 


3 


1  -  v 

1  +  V 


K 


K  ^ 

v 


T, 


(B. la) 


a 

n 


3v 

1  +  v 


K  »: 


K  a  T. 

v 


(B. lb) 
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The  ususl  ecttvsnSifis  of  2  2tr*s®  h»<n«  no*i tiva  tn  tension  has  been  used . 
The  various  quantities  appear ingare  Poisson's  ratio  v,  the  bulk  modulus 
K,  and  the  volume  coefficient  of  thermal  expansion  &  *  When  v  =  1/2 %  the 
coefficients  of  the  strain  are  both  equal  to  1  and  there  is  no  difference 
in  the  stresses* 

The  pressure  is  now  defined  by 

c  +  2^  V  -  V 

P  =  -  — 5 - -  =  -  K  «■:  +  R  »v  T  -  K  — - — -  +  K  <*v  T.  (B.2) 

o 

The  thermal  expansion  term  of  this  expression  is  associated  with  the  inter¬ 
nal  energy  term  of  Eq.  (3.35)  and  hence, 

K^T-^  G(r?)E,  (B.3) 

Eliminating  the  strain  between  Eqs.  (B.la)  and  B.2)  and  making  the  sub¬ 
stitution  of  Eq,  (B.3)  yields 

°  =*2tt^p  +  2tt tg^e  <»»•*> 

where  O'  is  still  positive  in  tension.  For  v  «=  1/2,  <J  »  -  P,  as  it  should. 

Eq.  (B.4)  is  suggested  as  the  proper  relation  between  the  pressure  of  the 

equation  of  state  formulation  and  the  stress  which  determines  material 
motions  whenever  this  difference  is  taken  Into  account. 
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APPENDIX  C 

EQUATION  OF  STATE  SUBROUTINE  LISTING 


Section  4  summarizes  the  equations  incorporated  into  the  equation  of  state 
subroutine,  enumerates  the  calculations  necessary  to  determine  certain 
parameters  and  to  define  boundaries  in  the  phase  plane,  and  indicates  a 
possible  checking  procedure.  It  also  lists  the  required  input  data  with 
units  which  must  be  made  available  to  the  program. 

Following  is  a  listing  of  the  program  which  determines  the  pressure  in  kbar 
for  a  given  volume  in  cc/gm  and  energy  in  cal/gm.  The  program  is  written 
in  Fortran  IV  for  a  time-share  system.  It  may  be  easily  adapted  to  Fortran 
IV  for  other  computers. 
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I*M 


-  .  »G  » I N 1 
1040  REAL  1NIT,K0,K1 ,K2,La,MU 
KWO  FF  C  MU)  =  MU*MU*(D+S»MU) 

- ffirfTTIl  "=  “5* 

70  FFP("!U)  =  (  D+  D+3  , *  S*  MU ) *  MU 

BO  GQP (ETA)  =  G 


nn  1,300, 

1100C  INITIAL  PNCE-CNLY  CALCULATIONS 
1110  300  ACC  =  IE-6 


1  =  1,93 
R.'  :  .0320597 
VO  =  1 .  /i;o 


1150  IF ( DSn)  1,1,2 

1160  l  FSO  :  WSf1*lE3/WT 
1170  SPSO  :  ESO 


3051  TP  J80 
GO  TF  3000 
8051VS0  s  VS 


USO  =  l./VS 
or  TO  3 
03  0  =  1)8  M 


0  :  1,/DSO 
3  ELO  =  E30+1  E3  *  (HLM-K5M) /WJ 
IF  ( 1)1. M)  A,  A  ,5 


0  =  .935* 

OF  TO  G 
5  DLO  =  jLP 


0  =  1./01.,) 
T"1  :  TMK/TCK 
Tl:  :  H;K/TCK 


c  c 

UL  s 
_ . 'ey  s 

l E3*IHLB-HLM> /(T3K-TMK) 

.  l£3»vnV*B-HVTiJ  /(73X-7MR)  *  . 

4 

_ 4-  .  - 

1420 

.--4-416-- 

1  440 

TJUTC 

AZ.X. 
4  1  = 

“TTV=CT -  - ' - 

-LLHsIILICJlI  BO/lRi*.l£lU 

PLTC/F  1 

- - - - - 

1  450 

ALPHA  z  A 2  +AI 

146  3 

y  z 

36./Tj  +  f..»T:.**6-42, 

1470 

AR  = 

(A2/TB+A1+.31425*X7 /<  l.+.J338*X> 

rvn~ 

75  : 

~T~.  R  7  .'7  r +TTT^Ti C\T;  -  771 ) 

1  490 

R(v  i 

3=1. /3 . 

1500 

RL  = 

1  .-‘■Cl  *  (  1  , -7f*.)  »*P?P72+3l  *(  1 

1510 

DC  s 

TOT/ftL 

1520 

VC  2 

1  ,/uc 

1530 

PC  = 

2C*R2*TCK*DC/VT 

1540 

B 1  : 

3  . 

1550 

7 

no  : 

31 

15  CO 

b\  - 

3  .«(  1  1  . )-(  1  . /7.C+3.  +  1  « /faO  >  /BO 

1570 

IF(ABS( CE1-B0) /B1 )-ACC)  8,3,7 

I5R0 

3 

B  = 

<<3.»BI-6.)*B1-1.)/CBI*<3.*BI-1 .)> 

1590 

BP  = 

(B 1 -3 , > /C3  ,*B  1*1.) 

160  0 

1F(K0)  10,9,10 

1610 

S) 

KO  r 

[J1 

1520 

10 

K  1  : 

B 1  -K  0 

153  0 

K2  = 

( 1 .+K 1+51 -ALPHA) /2, 

1540 

EPS  1 

=  ZC*R1*TCK/WT 

mo 

EPS2 

s  TCK*<CV-R 1 ) /WT 

1650 

EO  = 

( 1 ES*h  VB-C V* T3K ) /WT 

1670 

T  = 

TM 

1630 

ASK  I 

G  N  90  5  1  Tf  JVO 

1690 

GP  TP  9000 

1700 

905! 

1  Cl  = 

EO+ELO-EL 

1  71 0 

EVQ 

2  EV+El-EO 

1720 

PVO 

=  P~V 

173  0 

DVO 

:  l.'OPV 

1  740 

WO 

=  l./DVO 

uw 

££>  = 

El 

1760 

EC  = 

E0+EPS2-EPS1 *(KO+K 1+K1-K2) 

1770 

DLTEO  =  ELO-ESO 

1780 

EE  : 

A  •  *  ESQ  *DL  TEO 

1790 

DLTVO  =  VLO-VSO 

1  ROD 

f)  c  V 

:  1H  TEO/'I.TVO 

mo 

EQ  VO 

:  .04  l?.93*DLTEO/oLTVO 

1 920 

EES 

=  VO*  DEOV 

1330 

EET 

=  £0 VO*  EES /ESO 

1^40 

Z'u'J 

-  ES0*DLTVJ 

1W0 

CCz 

(<C1 /01 )*-3) /27 

1^60 

C  0  J 

:  1  £3  »HSM/(  WV  * (  TMK- J. 98, ) ) 

1370 

CEO 

sCL/wT 

13H0 

C3= 

(CL0+CS01/2, 

1390 

DLTCOsCLO-CSO 

1900  DCT  =  1 LTCO*  Tt'K 
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i9io  k  c  1  =bc  r /< . 

IDl’O  CE-TsC.  *  TVK 

1930  l  '• : El.j-C  50*  TMK 

»«»/./*  :  ”  ""  . . . 

t  -/  -» v  v  w  a  i*m 

1950  Yl=2.*Ct>T 

.ii£0  Y3.SYl*DCr _ _ 

1970  IN  17=1 

-f.P  Pr  i  *J  I  T I  a  .  C  A  LCtILA  » 1 1  *JS 
19*00  Giv-V;  r'I5;  -EG  I  PG  if  PHASE  UIAGPAK 

Sooo  ill  IF(E-ESO)  i5o,i5o,ii£ 

2010  112IF<V-VLO>130,!34,  139 

2020  130  Gg  Tg  4000  _ 

2030  4051  I FC V-VS) I  50 , 150, 131 
?040  131  lF(E-ELO)  132, 13?., 134 
2050  132  rf-,A>:  :  <-> E.30  )  /  (  EL3-ES0  ) 

2(te0  VKAXfc  =  FMAX*VLO+C1.-FMAX)*VSO 

2070  IF( VMAXE-V) 984 ,934,133 

2060  964  P  =  0, 

2090  GP  TP  9998  ‘ 

2100  133  GP  TP  4500 

2110  455121  TP  ISO 

Z120  134  5flTiT50o0 - 

2130  5051  IF(V-VLM) 133, 160,135 
2140  135  IF(E-EC)  136,  137,138 

50  136  GP  T?  5500 
21  SO  555  1  GP  TP  170 
2170  13  7  VLB  =  VC 

~2T30 - Pir"=  pc*n.-pvo> - 

2190  60  TP  170 

2&Q0  133  02  T0  6000 
'*5Ho  50  51  GO  TP  170 
2220  139  JF(E- EL 0)132, 132, 140 
2230  1 40  IFC V-VC) 141,144,146 
224b  i4i  Of  T0  65 oo 
2250  6551  IF (E-ELS) 1 42, 180,134 

2260  142  GO  TC  7000 _ 

2270  7051  GO  TP  130 
2230  144  SLB  =  EC 

2290  PV  =  1  . 

TSoo  TT~  i. 

2310  IF(E-EC)142, 180,190 

2320  146  IF(V-VVO) 147, 190,190 _ 

2330  147  GP  TP  7500 

2340  7551  If  <  E-EV)  142,180,190 

2350C  END  gf  PhASE  LOCATION _ 

2360C  SOLID  PHASE  P(V,E) 

2370  150  ETA  =  VO/V 

2330  m  =  ETA-1, _ 

ii'90  P . = . C *  MU+  FF  <  ;iU )  +G Ci  (  E T  A  )  *  E 

24QJ  GP  TP  9908 

2410C  Jl'LIJ-LICiUIL  MIXED  PHASE  P(V,E) 
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23150  WTf  =  "E0'vo*TTF- f . ) 

2430  GF  TO  9998 
9&&JJAUIP  PHASE  PCV.E), 
2350  1  70  RM=  I  .  /VLM 

2460  R8:|./VLB 


2460 

— .Q. 

2400 

2600 


2660  1 71  F =0 


Zr :  =  <  R i> P  L M-K  M*>'L3 )  /  (K  M -R  5  ) 

PI  =Z  1  /V-7  2 

16  =A  LOG  (PLM/PLB)  /AL0G(RM/RB) 


Z  4  =  ( A  |_f  0  (  RB  )  *  A  L  f’G  (  PLS  )  -  A  L  ’G  ( lr  ‘1 )  *A  |_  i’G  ( PLB )  )/AL2G<nM/Rb ) 
Al.F2  =  Z3*AL:’0(  |  ,/0)-7 

Ft  (Pl.M/(Rf1-LLJ)-Z3*F!.M/kM)/(7  I  -Z3*PLM/h  *1) 


IF(F) 171,174, 172 


GO  T0  174 


,-r  )* 


a&ymm  mtmjzxii 


( F+A  L0G( P 1 
2610  GF  T?  9993 

262 0C  LXQUIb-VAPCK  MIXEO  PHASE  P(V.£) 


ASSIGN  855  1  T0  035 

GB  TB  8600 


2690  8561  P  =  POPG 
2700  GF  TO  9 9 9 3 

.-’.71  00  PUGIN  OliUliFUTI  "EG 


1  1 

2730  4000  Y2=E-EB5 

2)40  TF=(Y2+SQRT<Y2*Y2«Y3)  )./Y  1 


50  ETA:  1 , 

760  HP  SO  :  E 

770  ASSIGN  SO 5?,  Tv"  J80 


ki;iu 


2790  3052  G0  TB  4051 

2800  C  GIVEN  V,E  IN  S-L  REG I  BN)  FIND  T(V,E) 


Kllll 


2850  4501  TFO: TF 
2860  ETAs VO/VS 


:  =  t  I  A  -  I  . 

0  ES?=C[*T-h:)C7/(TF«TF) 

0  ETA? : -ESP*  ETA  « E  f A /EES 


0V0*(TF-1.  )-C*MU-FF(riU)-GG(ETA)*ES 
HP:EOVO-(C+FFP( MU>+GGP(ETA)*ES)*ETAP-B3(ETA >*ESP 

. . _ _ _ : _ _ _ _  _ _  . 


.  030 
.  ’#*)U 

.TJVJ 

£9: trS+CB  r«'Ilr+HiJC  T/Tf 
vS:  <  ES-EPfi)  /vt.  jv 

I  f  C  A 1 [><  C  T  F- 1 r  0 )  /T  F)  -A  CC )  4500 . 45  JO  .  450  1 

2970C 

as’.o 

GIVEN  EL0«E:  FIND  VLM( E) ,PLM( E> 

5000  Y2sE-EBl. 

<1 9  JJ 

T  r  :  ( Y 2  +;.Vh  f  ( V2*Y  2 -V  3~> )  7?  1  * 

vj  J  V  C/ 

EH:  ?.:*S+Cb7*7f  +Kl)CT /Tr 

30  1 J 

ETA:  |  . 

3030 

ASSIGN  3033  TP  J30 

3040 

GE  T0  3000 

3050 

30*53  vimvs+cE-fes) TUTiTv 

3060 

FLHspS 

30  70 

G;'  Tr  505  1 

ibsoc 

Given  elo«e«ec*  find  vie< E) .FlSTZ) 

3090 

5300  Tl  s  TM 

3100 

TU  s  1. 

3110 

5501  T  =  (Tl+fl'j/2. 

3120 

ASSIGN  9052  T C  JO J 

3130 

ir  Tr  90U0 

3140 

9052  I FC AUS< ( E-EL) /£) -A CC ) 5505 ,5505 , 5502 

3150 

5502  IF(E-F.L)  55  03, 5505*5504 

3160 

5503  TU  :  T 

31  70 

or  T£  5501  ~ 

3  1  30 

5504  TL  :  T 

3 1  00 

g;  t.-  5501 

3&00 

5*05  VLB  =  VC/RL 

3210 

PL?>  =  P  C*(PV-PV0) 

3220 

G3  TE  5551 

323  OC 

GIVEN  EC«E|  FIND  VL;'.  :  VC  ,PLE (  E) 

3240 

6000  VL  ■  :  VC 

3250 

R  V  :  1. 

3260 

ASSIGN  3552  TP  J85 

' 

5270 

GE  T0  3500 

3?SQ 

8552  PLB  =  PC+PG 

3290 

Gr  T£  6051 

3300C 

:  GIVEN  VLQ*V«VCS  FI  NO  ElbCV) 

3310 

6500  KL  :  VC/V 

5320  XI  =  (  l  ,-RL)  / D 1 

3330  X  =  SQRT(XI*Xl/4.+CC) 

3340 _ X2  r  <X-X  1  /2  ♦  )**R00T3 

3350  X3  =  “(y+Xl/2.)»*R0«JT3 

3360  Y  =  X2+X3 

3370  T  :  1 .-Y +*3 


3  4*0  I  F £  ;-Vf 1AXE)  7001  ,7010.7010 

3400  7010  :  0. 


T5^5 - K*"T<rSF5?3 

3510  7001  T  :  (Tl.+  TU)  /2 . 

3520 _ ASSIGN  9054  Tfr  J90 _ _ 

PTiO  GO  T?  0300 

3540  0054  r  -  ( 1  , /R -  1 ,  /F.  L) /< 1  .  /F  V-  1 .  /R!  ) 

3553 _ EV  :  F«  EV+  ( ]  . -K  )  *  EL _ 

'3560  iFtAfiSl  < E-ED/E) -ACC >  7005  ,7005,7002 

3570  7032  IF< E-ET) 7003 ,7005,7004 

3580  7003  TU  =  T _ 

3r)!/0  j:  7  70  jl 

i'  OJ  7034  TL  r  T 

3r>  1  J _ 0  TO  700  1 _ 

5SfcO  7005  GO  TO  ToTT 

3630C  GIVEN  VC<V<VVO{  FINO  EV(V) 

3640  7500  RV  =  VC/V _ 

3450  X  I  :  FiV/<Z'>U  (E--fcP*hV)  *8V)  )-K2*RV«*3 

366  J  Xl  :  -K  0*r<  V*  ii  V 

3.i  70  X  .4  s  ( K  2  ■  V  “  K  1  )  *  n  V  *  K  V  _ _ _ 

3G80  Tal."  ’  '  '  ~~  '  ‘  . 

3600  PV  =  EXP<A2*Ct.-l,/T)-fAl*ALeGCT>) 

3700  7501  PVT  =  PV 
TfTo - Tim 

372  0  P'3  ;  X  1  * T+X2.+X 3 /  T 

3750  PVP  =  F V*  (A2/T+A  1  )/T 

T743  FSF  S  X 1“X3 /(T*T) 

3T50  T=TO+<PG-PV) /(PVP-PGP) 

3760 _ PV  s  EXPCA2*(1 .-1 ,/T)+A 1 «ALEG( T) ) _ 

TTT5  fF  ( A  n  S  (  C  P  V  -  P  V  T )  /P  V )  -  A  C  C  )~1 5  0  2 , 7  5  0  2  ,  T 5  0  1 

37R0  7502  EV  =  E0+EP82*T-F.PS l  ♦  C  (KO+2  .*K 1 /T)  -K 2*RV/T)  *r.V 
3  740 _ G£  TT  75  5  1 _ 

3^006  GlVfiN  SSO^EPSOt  FIND  VSCEPSO) 

38lO  8000  PS=E0V0*(TF-1 ,) 

3820  8001  ETAO  ■  ETA _ 

“3F30  '1U  :  if  A  -1  . 
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